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16 Abstract This report presents the resalts of the third and fourth year of work on a four-year 

research program which is part of the NASA HOST Program. The program goals were (1) to 
develop and validate unified constitutive models for isotropic materials, and (21 to 
demonstrate their usefulness for structural analyses of hot section components of gas 
turbine enaines. The unified models selected for development and evaluation were those 
of Sodner-Partom and of walker. 

The unified approach for elastic-viscoplastic constitutive equations has been shown 
to be a viable method for representing and predicting material response characteristics 
in the range where strain rate and temperature dependent inelastic deformations are ex- 
perienced. " This conclusion has been reached by extensive comparison of model calcula- 
tions against the experimental results of a sophisticated test program on two high tem- 
perature Ni-base alloys, B1900+Hf and Mar-M247, over a wide temperature range for a 
variety of deformation ana thermal Histories including uniaxial, multiaxial, and thermo- 
mechanical loading paths. Particular material effects which have been examined in this 
program include strain-rate and temperature dependence of inelastic flow, isotropic and 
directional hardening, static thermal recovery, multiaxial cyclic hardening, thermal 
history effects, strain aging effects, and the effects of thin coatings and grain size 
on mechanical response. Procedures for evaluating the material constants have also been 
developed and verified. 

The applicability of the Sodner-Partom and the Walker models for structural applica- 
tions has been demonstrated by implementing these models into the MARC finite element 
code and by performing a number of analyses including thermomechanical histories on com- 
ponents of the hot section of gas turbine engines and benchmark notch tensile specimens. 

The results of this four-year program have Deen documented in four annual reports. 
The results of the base program are summarized in this report. The tasks covered in 

this report include: (1) development of material test procedures, (2) tnermal history 

effects, and (3) verification of the constitutive model for an alternate material. 
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PREFACE 


This final report covers the effort of the optional program on the 
NASA HOST Program entitled "Constitutive Modeling for Isotropic Materials" 
conducted under Contract NAS3-23925. The effort of the base program is 
documented in the Second Annual Status Report. The NASA program manager for 
this project was Mr. Albert Kaufman who retired during the latter phase of 
this program. Dr. Robert Thompson has been the NASA Program Manager since 
Mr. Kaufman's retirement. The program manager at Southwest Research Institute 
is Dr. Ulric S. Lindholm. Contributors to this report are Drs. Kwai S. Chan, 
Richard A. Page and Mr. Andrew Nagy of SwRI, and Professor Sol R. Bodner of 
Technion, Israel, a consultant to this project. 
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1.0 INTRODUCTION 


The objective of the present program is to evaluate and further 
develop existing unified constitutive models for use in finite-element struc- 
tural analysis of turbine engine hot section components. This effort leads to 
a different approach for nonlinear finite-element computer codes which have 
heretofore been based on classical inelastic methods. The unified constitu- 
tive theory does not require many of the simplifying assumptions of classical 
elastic-plastic theory and is characterized by the use of kinetic equations 
and internal variables with appropriate evolution equations for treating all 
aspects of inelastic deformation including plasticity, creep, and stress 
relaxation. Specific model development has been directed toward isotropic, 
cast nickel-base alloys used for air-cooled turbine blades and vanes. The 
results of the present program have shown that this approach is particularly 
suited for determining the cyclic behavior of superalloy type blade and vane 
materials under cyclic loading conditions and high temperature environments, 
and is entirely compatible with three dimensional inelastic finite-element 
formulations. More efficient and accurate inelastic analysis of hot section 
components--turbine blades, turbine vanes, combustor liners and seals-- 
fabricated from "age-hardenable" isotropic superalloy materials has been 
realized as the result of these developments. 

The program was conducted in two phases: a basic program (Tasks A 
through I) and an optional follow-on program (Tasks J through M) . In the 
basic program which was completed in 1985 [1,2], two unified constitutive 
models, i.e., those of Bodner-Partom [3,4] and Walker [5], underwent further 
development for the prediction of the structural response of isotropic 
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materials for temperatures and strain ranges characteristic of cooled turbine 
vanes in advanced gas turbine engines (Task A). A data base of uniaxial and 
multiaxial material response characteristics required for the constitutive 
model development and evaluation was obtained for the base material, B1900+Hf 
(Tasks C and E). The constitutive models were then incorporated into the MARC 
finite-element computer code (Task D). The capability of the analytical 
method to predict the structural response was evaluated for multiaxial stress 
states (Task E) and nonisothermal conditions by conducting thermomechanical 
loading and benchmark notch verification experiments and analysis (Task F). 
As a final evaluation of the analytical model, a structural analysis was per- 
formed for a hot section component fabricated from the base material for 
simulated engine operating conditions (Task G). 

In the optional program, material property test procedures were 
developed further to obtain the material model constants from a minimum number 
of conventional tests (Task J). Further evaluation of the model was 
undertaken to consider thermal history effects (Task K). Finally, the 
developed constitutive model was applied and evaluated for an alternate 
material, Mar-M247 (Task L). 

The work under this program has been conducted as a joint effort be- 
tween Pratt & Whitney Aircraft (PWA) and Southwest Research Institute (SwRI) 
with technical assistance from Prof. Sol R. Bodner and Dr. Kevin Walker in the 
area of constitutive model development. The work and data base generated has 
been coordinated closely with another NASA funded HOST program at PWA (NAS3- 
23288) to develop advanced life prediction techniques for isotropic superalloy 
vane and blade materials. 
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The results of the present program have been reported annually in a 
series of NASA reports [1,2,6]. Results of the base program are summarized in 
the Second Annual Report, while this report documents those of the optional 
program. Some of the work in the optional program has been reported earlier 
in the Third Annual Report; part of these will be repeated in this report for 
the sake of clarity and completeness. While the base program demonstrated 
that both the Bodner-Partom and the Walker constitutive models were able to 
describe the extensive material data base with acceptable accuracy and were 
demonstrated to be efficiently utilized in the MARC finite element structural 
analysis code, the effort in the optional program has dominantly focused on 
the Bodner-Partom model for further study. This is not based on the advocacy 
of one model over the other but on the limitations of time and resources. 
Both models have certain advantages as well as limitations. Both will find 
further applications. 
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2.0 DEVELOPMENT OF MATERIAL TEST PROCEDURES (TASK J) 


2.1 Determination of Material Constants 

The procedure initiated in the base program for determining the 
Bodner-Partom model constants was developed further in Task J of the optional 
program. As a consequence of these efforts, the amount of required testing 
has been minimized. With the present methods, it is possible to determine the 
material constants of the Bodner-Partom model from conventional tensile data 
with supplemental creep data if slow strain-rate tensile data is not avail- 
able. Since a detailed presentation of the method for constants evaluation 
has been presented in the third annual report [6], only a brief outline is 
presented here. Procedures for determinating the material constants in the 
Walker model were established in the base program. Summaries of the Walker 
model, the procedures for determinating model constants, and material 
constants for the BI900+Hf alloy were presented in [2]. 

A summary of the Bodner-Partom model is given in Table 2.1, while 
the material constants in the model are listed in Table 2.2. As illustrated 
in Figure 2.1, the proposed method relies on the use of work hardening data in 
the form of a y versus a plot to obtain the hardening constants, where y is 
the derivative of the stress with respect to plastic work (y = do/dW p ). An 
example of a y-o plot for B1900+Hf is shown in Figure 2.2. The y-a plot pro- 
vides values for and m 2 (assumed to be constants), and also the value of 
the saturation stress, , at the test strain rate. The dependence of on 
strain rate is used for determining the value of n. Figure 2.3. It should be 
noted that the required ct $ data can be generated by either performing tensile 
tests on multiple specimens at various strain rates or on a single specimen 
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TABLE 2.1 


A SUMMARY OF BODNER-PARTOM MODEL 


Flow Law: 


. e 

+ e P 

= e iJ 

e ij 

3 x S i 

j ; d kk = 1 

' a i j 

‘ I °'ij°kk 


Mj 


Kinetic Equation: 

D ? 

D 2 * D Q exp 

with Z = Z* + Z D 

. P 1 . p . p 
°2 3 2 e 1j e ij 

J 2 = I S 1j S ij 
\ z * d 2 p /j 2 



Evolution Equations of Internal Variables: 
a. Isotropic Hardening 


Z 1 = m, aZ 3 - Z 1 j W p - A 1 Z 1 2 i | 


"1 ^1 3 "3 


n 


where d = m 2 C a-, sine - ct)W p 
0 = COS 


1 (V„V„) or 0 = cos' 1 (u.j u.j) 


ij iJ' 


V 1j * 6 ij^ 3 kl s kl^ » v ij 3 ^i j/'®kl®kl 1 
u ij 3 ff ij /(o kl a kl )1/2 . “,-j 3 ^ij/^kl d kl )1/2 

,1 - ■ . , . P 


with Z 1 (o) = Z 0 ; fl p = ; VJ p (o) = 0; a(0) * 0 

b. Directional Hardening 

®1 j 3 m 2 (Z 3 U ij - V *p • Vl 


(BklSkl) 


1/2 


V U 


W 11 


h Z D - 0 i ,-u . . ; Z D (o) = 0 , B,,(o) = 0 

I U 1 j 1 J 


Material Constants: O 0 , Z Q , Z] , Z 2 , Z 3 » . m 2> G 1 

A-| , A 2 , r-i, r 2 , n, ar.d elastic constants 


In most cases can set: r-j = rg, A-] = A 2 , Z Q = Z 
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TABLE 2.2 


MATERIAL CONSTANTS IN BODNER-PARTOM CONSTITUTIVE MODEL 


Constants 


D 

Z 


o 

o 



Descriptions 

Limiting shear strain rate, sec" 1 (assumed value) 

Initial value of the isotropic hardening 
variable, MPa, (T) 

Limiting (maximum) value of Z 1 , MPa 

Fully recovered (minimum) value of Z 1 , MPa, (T) 

Limiting (maximum) value of 2 ®, MPa 

Hardening rate coefficient of Z 1 , MPa" 1 
Hardening rate coefficient of Z D , MPa" 1 

Kinetic parameter, (T) 

Recovery coefficient for Z 1 , sec" 1 , (T) 

Recovery coefficient for Z*^, sec” 1 , (T) 

Recovery exponent for Z 1 , (T) 

Recovery exponent for Z* 1 , (T) 

Limiting coefficient for non-proportional 
hardening effect 


Note: Symbol (T) indicates general temperature dependence 
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Tensile Data 


fit a and to a pclynominal and 

obtain v = - -AS (tensile data) 

a aen 


obtain and from the slopes 
of y vs. a (tensile data); set D Q 


obtain n from the slope of saturation 
stress vs. strain rate (tensile data) 


calculate the sum of Z. 

and Z^ from saturation 

stress, n and i p 

(tensile data) 


obtain Zg from .2* offset yield 
stress; set = Z Q 


calculate from yield and saturation 
stresses (tensile data), obtain Z3 


calculate A and r from slow rate tensile data or 
steady state creep rate f rom q vs. (^ c “^ 0 ) plot 


Tensi le Data and/or 
Creep Data 


FIGURE 2.1 A FLOW CHART SHOWING A SYSTEMATIC PROCEDURE 
FOR EVALUATING THE MATERIAL CONSTANTS IN 
THE BODNER-PARTOM THEORY 
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FIGURE, 2.2 THE ISOTROPIC AND DIRECTIONAL HARDENING COMPONENTS OF B1900+Hf 
IN A Y-o- PLOT 
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FIGURE 2.3 DETERMINATION OF THE VALUE OF n BASED ON 
SEVERAL STRAIN RATES 
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involving strain rate jumps. In the latter case, the strain-rate jump should 
be made after saturation has occurred. Figure 2.4 indicates that the satura- 
tion stresses at e - 1 x 10"^ sec - ^ obtained based on these two approaches 
show little difference. Once n is obtained, the values for Z Q , and Z 3 can 
be computed from the intercepts in Figure 2.2 using the 0.2% offset yield 
stress (or the proportional limit), the saturation stress, and an assumed 
value of D q . The procedures outlined so far are applicable for temperatures 
at which static thermal recovery is unimportant. For elevated temperatures at 
which static thermal recovery is important, the thermal recovery constants 
need to be evaluated. This is accomplished through the use of a plot of q 
vs. ° c -° 0> where 

q = (L,|m I (Z r Z 0 ) + «I 2 (Z 0 +Z 3 )| - m 2 o c )„ c ^ (2.1) 

» 0 * [2tn(2D 0 /z3 e p )|-‘ in Z 0 

and o c and are the stress and plastic strain rate under steady-state condi- 
tions. A plot of q vs a c -a Q for B1900+Hf is shown in Figure 2.5. As dis- 
cussed in [6,7], the slope of such a plot is the thermal recovery exponent r, 
and the values of q and ° c - o 0 can be used to compute the thermal recovery co- 
efficient A. This procedure is based upon assuming that Aj=A 2 and r 1 =r 2 =r. 

The proposed procedures have been applied to the B1900+Hf and Mar- 
M247 materials in this program with good results. They have also been applied 
to HY-80, A533B steels, and Astroloy 700 in other programs at SwRI. In all 
cases, very satisfactory results are obtained. 

A limited amount of effort has been devoted to examining the possi- 
ble effects of grain size and surface coating on the model constants through 
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FIGURE 2.4 COMPARISON OF MONOTONIC TENSILE AND JUMP TEST RESULTS OF B1900+Hf 




MPa/sec 



5 10 50 100 500 

V V Mpa 


FIGURE 2.5 DETERMINATION OF THE VALUE OF r BASED ON 
B1900+Hf ES ° F THE q VS °c"°o CURVES F0R 




their influence on the constitutive behavior of the alloys. Figures 2 . 6 - 2. 8 
show that for the B1900+Hf material investigated, grain size has only a small 
effect on the stress-strain behavior under either tensile or cyclic loading. 
The effect of a thermal barrier surface coating on the deformation behavior of 
B1900+Hf was investigated by performing a TMF test and those results will be 
presented in the next section. It was also found that this coating also had a 
negligible effect on the model constants. 

2.2 Strain-Aging Effects 

Step strain-rate tensile tests in the optional program indicated 
that negative strain rate sensitivity, which is a consequence of strain aging 
effects, is present in B1900+Hf at a temperature range bounded above by about 
760°C, and also exists in Mar-M247 in the temperature range of 23-648°C. 
Under this circumstance, certain predictions of the Bodner-Partom model do not 
correlate well with the experimental data because strain aging effects are not 
included in the model. In addition, it has been suggested that strain aging 
can lead to thermal -hi story effects which can potentially invalidate isother- 
mally based constitutive models subjected to nonisothermal loading histories 
[8,10]. As a result, a small effort has been attempted to incorporate the 
strain aging effect into the Bodner-Partom model in order to improve the 
accuracy of the model in the regime of influence. A strain aging term is 
present in the Walker model [2]; thus, modification of this model is 
unnecessary. 

Strain aging effects are primarily due to the interaction of dis- 
locations and solute atoms with moving dislocations [11-13]. Common observa- 
tions associated with strain aging are: (1) serrated flow, (2) negative 
strain rate sensitivity, and (3) a plateau in the temperature dependence of 
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FIGURE 2.6 COMPARISON OF JUMP TEST RESULTS FOR B1 900+Hf 
SPECIMENS AT 648°C AND WITH TWO DIFFERENT 
GRAIN SIZES. 
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FIGURE 2.7 COMPARISON OF JUMP TEST RESULTS FOR B1 900+Hf 
SPECIMENS AT 982 °C AND WITH TWO DIFFERENT 
GRAIN SIZES. 
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the yield stress. In a review paper, Schmidt and Miller [12] indicated that 
solute enhanced strain hardening can be due to a number of processes which 
include precipitation and the solute affecting the rate of dislocation 
multiplication, annihilation, or the recovery of dislocation structures. As 
indicated earlier, strain aging and precipitation have also been linked to 
thermal history effects [8-10], 

Schmidt and Miller [14] have incorporated a strain aging term into 
the unified constitutive equations of Miller [15] by introducing a function 
which is added to the isotropic hardening variables. The functional form pro- 
posed by Schmidt and Miller to represent solute strengthening is: 

F(solute hardening) = a Q exp 1-aj log(e/e Q ) 2 ] (2.2) 

where a Q , and a^ are constants and e and e Q are the applied and reference 
effective strain rates. It is noted that this procedure is an empirical 
method for representing the effect of strain aging and does not correspond 
directly with physical processes. This functional relationship was adopted in 
the Walker model for incorporating strain aging effects [2]. Walker, however, 
considered solute hardening to be directional in character and the strain 
aging effects were incorporated in the evolution equation for the equilibrium 
stress [2]. A possible method for modifying the Bodner-Partom model to in- 
clude strain aging effects would be to adopt Equation (2.2) into the evolution 

equation for the directional hardening variable, e.., so that 

J 

S i j = m^fz^ U ij ~ e ij (l -a )l - (thermal recovery term) (2.3) 

where a represents the strain aging term given by Equation (2.2). The 
constants a Q , a^, and e Q must be evaluated from tensile data in the strain 
aging regime. Figures 2.9(a) and (b) show the stress-strain results by 
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(a) 



(b) 


FIGURE 2.9 


10DEL REPRESENTATIONS OF B1900+Hf AT 649°C USING A 
STRAIN AGING TERM: (a) EFFECTS OF a p IN THE STRAIN 

\GING TERM ON THE STRESS-STRAIN CURVE, AND (b) STRAIN 
\GING EFFECT SIMULATED WITH THE BODNER-PARTOM MODEL 
JSING a = .5, aj = .1, and i Q = 1 x 10 sec 
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systematically varying the parameters a Q , aj, e Q , and strain rate using the 
material constants for B1900+Hf at 649°C. The results indicate that a 
strength peak can be produced at a strain rate of e = e so that the flow 
stress decreases when the strain rate is either increased or decreased, which 
simulates one of the characteristics of the strain aging effect. The 

magnitude of the flow stress enhancement depends on the value of the material 
constants a Q . Calculations of tensile curves with several strain rate 
changes, however, indicate that Equation (2.3) would not produce a strain-rate 
independent inelastic flow behavior. Although a negative strain-rate 
dependence is obtained, transient overshoots in stress, not observed in 
B1900+Hf, are often accompanied with the strain rate changes. Figure 2.10. 

Another characteristic of the strain aging phenomenon is the in- 
crease in strength with increasing temperature, resulting in a plateau or 
strength peak in the temperature dependence of yield strength. Such a 
strength peak exists at * 648-760°C in B1900+Hf and Mar-M247. However, the 
strength peak in this intermediate temperature range is generally attributed 
to the deformation characteristics of the Y ' phase and is not due to strain 
aging [16]. It should be noted that a strength peak is also observed in 
single crystal and polycrystalline N i 3 A 1 ( y' ) at the intermediate temperature 
range. The mechanisms of the strength peak in Y ' appears to be dislocation 
interactions resulting from cross-slip from ( 111 ) planes to ( 010 ) planes 
[16]. The latter type of dislocations are sessile and they tend to interact 
with subsequent dislocation movements on the (111) planes. The consequences 
are anomalous temperature-strength behavior, asymmetric flow, and unusual 
hardening behavior. Since B1900 + Hf and Mar-M247 contain approximately 60% 
(±5%) volume fraction of r, it is possible that the observed temperature- 
strength behavior and thermal-history effects in Ni-alloys are caused by N i 3 A 1 
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FIGURE 2.10 STRESS-STRAIN CURVE CALCULATED WITH THE B0DNER-PART0M MODEL 
WITH A STRAIN AGING TERM 





rather than strain aging* Because of the uncertainties of the underlying 
mechanisms, further attempts to include strain aging effects in the Bodner- 
Partom model were discontinued. 

Unusual cyclic hardening effects, to be presented in the next sec- 
tion, were observed in TMF cycling of B1900+Hf in the temperature range of 
538-760°C; transmission electron microscopy revealed that the anomalous 
hardening behavior originated from dislocation interactions in the Y ' phase. 

2.3 Nonproportional Hardening 

As discussed in [61, the Bodner-Partom model somewhat overpredicts 
the stress response of Mar-M247 subjected to 90° out-of-phase tension-torsion 
strain cycling. Similar results were also observed in the calculations for 
B1900+Hf [2]. In order to examine this effect, exercises were performed using 
the Bodner-Partom model and the constants for Mar-M247 at 982°C. The loading 
sequence was: (1) uniaxial strain cycling under an axial strain range of ±.456 

for 5 cycles, (2) 90° out-of-phase tension-torsion for two and one-half cycles 
at an effective strain range of ±.456, and (3) return to axial loading as 
described in Step (1). The calculated a vs. /3 t responses associated this 

sequence is shown in Figure 2.11(a). It is noted that cyclic saturation is 
almost attained after uniaxial cycling for 5 cycles. Additional cyclic 
hardening, however, is observed when the 90° out-of-phase cycle is applied. 
This nonproportional hardening saturated rapidly, and the o-/3x locus becomes 
stable in the second cycle. Detailed examination of the internal variables 
indicated that the calculated nonproportional hardening effect is the conse- 
quence of an increase of the directional hardening variable, Z D , when the 
loading condition is changed from uniaxial to out-of-phase loading. This 
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FIGURE 2.11 BODNER-PARTOM MODEL CALCULATION FOR MAR-M247 UNDER 

(1) UNIAXIAL CYCLING FOR 5 CYCLES, FOLLOWED BY 

(2) 90° OUT-OF-PHASE CYCLING FOR 2-1/2 CYCLES, AND 

(3) UNIAXIAL CYCLING FOR 5 CYCLES: (a) o-/3x PLOT, 
AND (b) TIME HISTORY OF Z D 









nonproportional hardening effect is transient and diminishes as the loading 
condition is returned to uniaxial loading. The corresponding history for the 
directional hardening variable, Z^, is shown in Figure 2.11(b). 

A possible explanation for this transient nonproportional hardening 
is that the conditions for cyclic saturation under uniaxial loading and non- 
proportional, multiaxial loading are different. For the latter case, it 
appears that 


Z = s. .u . . = constant 

1 J ' J 


(2.4) 


during cyclic saturation under a 90° out-of-phase strain cycle. As a result. 


Z = 6 . .0. . + 6. -U.. = 0 


(2.5) 


Z 


I 


= 0 


( 2 . 6 ) 


would be the conditions for cyclic saturation under nonproportional loading. 
On the other hand, the condition for cyclic saturation under uniaxial strain- 
ing is Z^ = 0, but Z ^ varies continually with cycling and its amplitude 
depends on the cycle strain range. 

This type of nonproportional hardening effect appears to be present 

in B1900+Hf and Mar-M247. The effect, however, is quite small and becomes 

saturated within one cycle. The overprediction of the axial stress by the 

Bodner-Partom model for 90° out-of-phase loading appears to be primarily due 

to overestimation of this transient, nonproportional hardening effect. 

Further examinations of the calculated results and Equation (2.5) revealed 

that the discrepancies can be completely eliminated by an adjustment of the 

thermal recovery constant A 2 in the expression fore.., as illustrated in 

J 
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Figure 2.12. Since the magnitude of directional hardening B-j is essentially 
constant under 90° out-of-phase loading, it would be sensitive to the thermal 
recovery constant at elevated temperature. It may therefore not be proper to 
set A 1= A 2 when dealing with multiaxial nonproportional loading at elevated 

temperature. 

2.4 Update on Finite-Element Analysis of the Benchm ark Notch Specimen 

Finite element analyses of the benchmark notch-round specimens were 
performed in the base program using the MARC Code [17] and the Walker consti- 
tutive model. The Bodner-Partom model was not used in these numerical appli- 
cations. In the last two years, there has been an effort at SwRI to develop 
an in-house elastic-viscoplastic finite-element code which is suitable for 
high strain rate (> 10 4 sec -1 ) structural applications. This finite-element 
code, which has been completed during the last year, uses the Bodner-Partom 
model with isotropic hardening. Because of the availability of this new 
finite-element code, an analysis of the benchmark notch-round specimen sub- 
jected to a monotonic load was performed. The analysis was intended as a 
check of the finite-element code and was carried out in another SwRI project 
at no cost to NASA. Since the results are of interest to this program, the 

benchmark notch calculations are reported here also. 

Figure 2.13 shows the finite-element mesh used in the benchmark 
notch calculations. It consists of 283 axis-symmetric elements with 1832 
active degrees of freedom. Since the benchmark notch test was conducted under 
displacement control, the load-line displacements were prescribed on the 
boundary elements. The finite-element calculations were performed on a 
VAX/780 computer, and the CPU time was approximately one and one-half hours. 
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TIDE <SFC> 
(b) 


FIGURE 2.12 BODNER-PARTOM MODEL CALCULATION FORMAR-M247 UNDER_THE 
SAME LOADING HISTORIES AS IN FIGURE 2.11; (a) 0 -/ 3 x 
PLOT, AND (b) TIME HISTORY OF Z D DURING CYCLING. Re- 
sults are obtained by increasing the value of A 2 by 
a factor of 10. 
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FIGURE 2.13 FINITE-ELEMENT MESH OF THE BENCHMARK 
NOTCH-ROUND SPECIMEN 
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Figure 2.14 shows the excellent agreement between experiment and finite- 
element calculation based on the Bodner-Partom model. Previously reported, 
the finite-element calculations based on the Walker model and MARC code 
overpredicted the applied load at a given diametrical displacement. These 
finite-element results, together with previous results on the turbine blade 
analysis, indicate that the unified constitutive models are viable for 
structural applications based on the finite-element method. 
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FIGURE 2.14 COMPARISON OF MODEL PREDICTIONS WITH EXPERIMENTAL 
RESULTS OF THE BENCHMARK NOTCHED SPECIMEN UNDER 
MONOTONIC LOADING 
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3.0 THERMAL HISTORY EFFECTS (TASK K) 


3.1 Thermal History Effects— Literature Survey 

A literature survey was conducted as part of Task K1 to assess the 
effects of thermal history on constitutive behavior and their implications 
with regard to constitutive modeling based on the unified approach. The 
literature survey indicated that the inelastic response of metals [8,9,18] in- 
cluding yield surface determinations [19], depends, in general, not only on 
the current temperature but also on the thermal history. This is analogous to 
the fact that inelastic deformation depends on the deformation history as well 
as on the current stress. In some instances, the hardening response of mate- 
rials is observed to differ under isothermal and nonisothermal conditions 
which has generally been attributed to either strain aging [8,9] or other 
microstructural effects [18]. 

Thermal history effects are usually modeled in the unified consti- 
tutive approach by including thermal terms in the evolution equations for the 
isotropic and directional hardening variables [10,20]. The general forms of 
the evolution equations for the isotropic hardening variables, K, and the 
directional hardening variable, a.., are as follows [20]: 


* = F(J 2 ,K,T) = h 1 (K,T)M 1 - r L ( T,K) 


(3.1) 


G. . 
1J 


- ^(o.j.n.j.T) h 2 (G. j ,T)|5|. j . - d(o 1 j,T)fi i j - r 2 ( n ij» T ) V ij ( 3 * 2 ) 


for isothermal loading, and 

* = h^K.T)^ - r^T.K) + e^K.TjT (3.3) 



h,(a. .,T)M. . - d(a. .,T)fi. . - r 0 (a. .,T)V. . 
2 V ij* ' ij v ij* ' ij 2 V ij* ’ ij 


(3.4) 


e 2 (G ij .,T)T W.j 
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where h^, r^, and 8^ represent, respectively, the hardening, static thermal 
recovery, and thermal history functions for K; d, r^, and 0 2 represent the 
hardening, dynamic recovery, static thermal recovery, and thermal history 
functions for n.^, respectively; — and N.j are the hardening measure; 
and are directional indices. Specific forms of ft.j, ^ j V-jj* and W^- are 
given in [20] . 

The appropriate forms of 9 1 and 0 2 are not well established at this 
time. A general approach for modeling thermal history effects is to express 
and 0 2 as functions of the internal variable and temperature [10,20]. New 
internal variables may also be introduced [18]. In a particular approach 
[2,5,21], 0 1 and 0 2 are taken as functions represented by the temperature 
derivatives of the evolution equations for isothermal conditions [2,5,21]. 


01 (T) = F(J 2 ,K,T), t 

(3.5) 

= G( fl i j’T, a i j) »t 

(3.6) 


which generally contain terms involving the temperature derivative of material 
constants. Temperature-dependent material constants in the kinetic equation 
are treated in a similar manner. Comparisons of thermomechanical constitutive 
data of B1900+Hf and calculations using the Bodner-Partom and Walker models at 
SwRI and PW&A indicated that this approach led to reasonably good agreement 
between theory and experiment [2] suggesting that thermal history effects 
associated with TMF cycling can be predicted based on isothermal data alone. 
No additional temperature rate (T) terms other than those represented by 
Eqs. (3.5) and (3.6) were needed. Whether or not this conclusion is valid for 
other materials or thermal paths remains an open question. 
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Recent thermomechanical studies [8,9] report response characteris- 
tics under nonisothermal conditions which are substantially different from 
those manifested under isothermal circumstances. While these studies clearly 
demonstrate the presence of strain aging and thermal history effects on the 
behavior of certain structural alloys, the results do not contradict the pos- 
sibility of predicting thermal history effects by an appropriate constitutive 
model which uses material constants based on isothermal data. A systematic 
evaluation of model prediction against experimental results obtained under a 
variety of thermomechanical paths is required to resolve this issue. 

Based on the results of this literature survey, SwRI recommended 
that efforts in Task K2 be directed to examining thermal history effects in 
B1900+Hf by thermomechanical cycling: (1) at a temperature range where strain 
aging effect is prominent, (2) at a temperature range where microstructure 
changes such as precipitation or resolution of precipitates could occur, (3) 
at various heating rates, (4) using different thermal paths, and (5) with 
elastic as well as inelastic paths. The experimental efforts have been 
coupled with parallel analytical efforts to determine whether or not thermal 
history effects can be predicted using the Bodner-Partom model in conjunction 
with material constants obtained solely from isothermal data. 

3 * 2 Thermal History E ffects: Numerical Integration 

The various unified constitutive models proposed in the literature 
contain material constants of which some are temperature dependent. These 
temperature-dependent constants can include the initial and, in certain 
models, the limiting values of the hardening variables, and also terms in the 
kinetic equation. Using these constitutive models for thermomechanical 
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loading requires a judicious means of tracing and updating the current values 
of the hardening variables and the other materials constants during tempera- 
ture changes. A theoretical framework for this purpose is presented in this 
section. The results of this investigation also indicate the general forms of 
evolution equations for the hardening variables when temperature rate terms 
are included. 

At a temperature T and a temperature rate T (T = dT/dt), all of the 
unified constitutive models can be represented in the following form: 

£ P " l (° * * » T ) (3-7) 

X = g (a, x , T) (3.8) 

o = 0 (e - e P ) + e(a , 0) T (3.9) 

30 

with 0 (a , 0) = [e - e P - a (T-T q ) ] 

- 9 ((T-y f? + °] ( 3 - 10 ) 

where e , e P , a , and x are the total strain, inelastic (plastic) strain, 
stress, and hardening variable vectors, respectively; D is the stiffness 
matrix and a is the coefficient of thermal expansion. The dot above the 
variables denotes the derivative of the variable with time. It should be 
noted that Equation (3.10) is obtained by differentiating Hooke's law with 
respect to time for thermomechanical loading conditions by including effects 
of temperature change on the elastic moduli and coefficient of thermal 
expansion. 
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Using Tanaka's approach [22], an implicit time operator is used for 
integrating the rate equations in the unified models. The incremental stress 
vector Ao, strain vector Ae, nonelastic strain vector Ae P , and structural 
variable vector ax, occurring in the time interval At = b t - a t can be defined 
as 


= At { ( 1-n ) + n b o} , 

(3.11) 

= At { ( 1— n ) a e + n b e} * 

(3.12) 

= At { ( 1-T1 ) a e P + n b e P } , 

(3.13) 

At{(l-n) a x + n b x} , 

(3.14) 


in which b ( ) and a ( ) represent the values of the parameters ( ) at 
t = b and t = a , respectively. The implicit time operator, n , has values 
ranging from 0 to unity, i.e., 


0 < n < 1 (3.15) 

It should be noted that the numerical scheme reduces to Euler's explicit 

method when n = 0, and to Euler's implicit method when n = 1. 

The values of e P and x , at time t = b can be approximated by their 
previous values and their gradients at t = a through Taylor's expansion: 

b e P = a f + a f Aa + 3 f AX + a f T AT (3.16) 

~,a ~,X ~ ~,T v ' 

b x = a g + a g -p Ae p + a g ax + a g at (3.17) 

~ ~ ~ s * ~ ~ 9 l 

Substituting Equations (3.7-3.10) and Equations (3.16-3.17) into Equation 
(3.11-3.14) and after rearranging terms leads to 
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(3.18) 
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-1 
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isothermal terms nonisothermal terms 


where 
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-nAt 


0 
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-nAt a f 




- a g 


l 9,X 
0 


0 

-D 


A V = 


At{ a f - n( a £ P - a f)} 


At{( a x - a g . D a E p ) - n ( a * - a 9> 


(3.19) 


(3.20) 


(3.21) 
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and 


r|At a f T f 

~ 5 ■ 



(3.22) 


It is important to note that Equations (3.19) are divided into isothermal and 
nonisothermal terms. The nonisothermal terms include the thermal load terms 
and the terms representing the temperature derivatives of the kinetic equation 
(f) and the evolution equations of the hardening variable (g). The 
temperature dependence of the material constants enters directly in the 
expression for a f^ T and a g^. The approach outlined above was applied to the 

Bodner-Partom model and used in all the numerical calculations shown in this 
program. 

For the Bodner-Partom model. 


f = x S. . 

ij 


g = e- • 


whose temperature derivative are as follow: 


~.T 


1 IJ. 

2 3J 


£n 


3J, 



(3.23) 


(3.24) 


(3.25) 

•p /dn \ 
i j l dT / 

(3.26) 
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Zm^/D,^ 
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T 



Vi 


i - 1. 




A t dT 



(3.27) 


q T = 2m 0 /D 0 J 0 u 
3 J 


2 /u 2°2 u ij V dT 



Vl 


n T* 


’ / 0 kfc B ka 


ill 


/6 ksi 6 ksi 



It is noted that this model contains relatively few temperature-dependent 
material constants (compared to other models) which simplifies the numerical 
work. From Equations (3.26-3.28) it is clear that the temperature derivatives 
of the material constants are needed in the nonisothermal load terms in 
Equations (3.19). The appropriate temperature derivatives of the Bodner- 
Partom model constants for B1900+Hf are obtained by fitting the temperature- 
dependent material constants, shown in Table 3.1, to appropriate functions and 
then taking derivatives of those functions with respect to temperature. Sum- 
maries of the expressions representing the temperature dependence of the model 
constants for B1900+Hf are given in Table 3.2. 


3 . 3 . Thermomechanical Constitutive Testing 

3.3.1 Experimental Set-up 

Thermomechanical response tests were conducted on an MTS hydraulic 
testing machine equipped with an inducting heating unit, a dynamic temperature 
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TABLE 3.2 


TEMPERATURE DEPENDENCE OF THE BODNER-PARTOM MODEL 
CONSTANTS FOR B1900+Hf 


n = 1.055 


T < 863. 4°C 


n = -1.0524 + 


2.394.74 


863. 4°C < T < 1093°C 


Z = Z, = 2700 MPa 

0 C- 


T < 838°C 


Z o = Z 2 = " 4950 + 8 (T+273) MPa 


838 °C < T < 1093°C 


Ai = A 2 = 1.319xl0 y exp 


30580 

(T+273) 


sec 


-1 


T < 1093°C 
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monitoring system, a thermal strain compensator, and a cooling system for the 
test specimen. A schematic of the TMF testing system is illustrated in 
Figure 3.1. The specimens could be subjected to controlled axial extension or 
compression at a prescribed temperature-time history. The dynamic temperature 
monitoring system, which consists of a number of thermocouples and isolation 
amplifiers, was constructed at SwRI and was used both for temperature control 
and for measuring the temperature profile along the gauge length of the test 
specimen. The thermal strain compensator, also built at SwRI, is an elec- 
tronic circuit device which reproduces the thermal strain of the test specimen 
at a given temperature. As indicated in Figure 3.1, the input to the thermal 
strain compensator is temperature (from the temperature controller), and the 
output is a voltage signal corresponding to the thermal strain of the test 
specimen at that temperature. In a TMF test, the thermal strain signal is 

added to the prescribed mechanical strain signal in the MTS machine 
controller. 

The thermomechanical test specimens were of tubular geometry, as 
shown in Figure 3.2(c), and were heated by induction heating and cooled both 
externally and internally by forced cold air (at RT). Internal cooling was 
accomplished by flowing room temperature air through the center of the tubular 
specimens and the grip system. The forced cold air flow was regulated with 
solenoid-operated flow distribution manifolds. The performance of the TMF 
testing system was evaluated by cycling from 538 to 982°C. Figure 3.3 shows 
good agreement between the prescribed temperature setting and the temperature 
on the specimen measured by the thermocouple located at the center of the 
specimen. The temperature distribution along the gauge length of the specimen 
was within 20 °C of the prescribed temperature. This temperature gradient was 
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FIGURE 3.1 A SCHEMATIC OF THE THERMOMECHANICAL TESTING SYSTEM 
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FIGURE 3.3 COMPARISON OF PRESCRIBED TEMPERATURE WITH SPECIMEN (MEASURED) TEMPERATURE 



less than 9°C during heating but a 20°C gradient occurred along the axial 
direction of the gauge length of the specimen during cooling. The lower part 
of the specimen was persistently hotter than the prescribed temperature 
because of the direction of internal cooling air flow. 

For calibrating the thermal strain compensator, the thermal strain 
values of B1900+Hf and Mar-M247 were measureu by heating the specimens, at 
zero load, from 38°C to 982°C and then cycling between 538°C to 982°C for ten 
cycles. Figure 3.4 shows that the thermal strain of B1900+Hf is a nonlinear 
function of temperature. Consistent thermal strains were obtained during 
thermal cycling indicating good reproducibility. The thermal strain data of 
Mar-M247 are identical to those of B1900+Hf. After calibration, the thermal 
strain compensator was evaluated by performing TMF strain cycling using the 
thermal strain compensator as the control input. Such a thermomechanical 
strain cycle should ideally lead to zero stress on the specimen. The maximum 

thermal stress observed on specimens under these check tests was less than 
14 MPa. 

Figure 3.5 shows a close-up view of a TMF specimen under test during 
TMF cycling at temperatures between 538 to 982°C. In addition to the 
induction coils. Figure 3.5 shows the small copper tubings used for external 
cooling and the externally mounted extensometer which is attached to the 
specimen through ceramic reach-rods. The reach-rods are placed on two small 
indents which are 25.4 mm apart and within the gauge length of the specimen. 
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FIGURE 3.4 THERMAL STRAIN OF B1900+Hf SPECIMEN HEATED FROM 38°C TO 538°C AND THEN 
CYCLED AT ZERO LOAD FROM 538°C TO 982°C FOR TEN CYCLES 
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FIGURE 3.5 A CLOSE-UP VIEW OF A THERMOMECHANICAL CYCLIC 
SPECIMEN UNDER CYCLING FROM 538°C TO 982°C. 

The specimen is heated by induction heating 
and cooled by feeding cold (RT) air internally 
through the tubular specimen and externally 
using small copper tubes. 
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3.3.2 


Hon- isothermal Loading Paths 


A number of experiments involving various thermomechanical loading 
paths were used to test the unified constitutive approach. The test matrix, 
shown in Table 3.3, consisted of simple as well as complicated paths of 
strain and temperature. Four types of nonisothermal tests were performed: 
(1) step temperature tensile tests; (2) step temperature creep tests, (3) step 
temperature cyclic tests; (4) thermomechanical fatigue (TMF) tests. The 
strain-time and thermal history for these tests are schematically illustrated 

in Figures 3. 6-3. 9. 

As shown in Table 3.3, six-step temperature tensile tests were per- 
formed. In these tests, the specimens were pulled at a constant strain rate 
up to a strain level of 1-2% at an initial temperature T^ With continued 
straining, the test temperature was then changed at a prescribed heating or 
cooling rate until it reached T 2 . After an additional 0.5%. strain at T 2 , the 
test temperature was changed again to T 3 in a similar manner. In some cases 

the specimen failed after only one temperature change. 

In the step temperature creep experiments, the specimens were first 
subjected to a constant stress oj. at temperature J l until a steady creep rate 
was reached. The test temperature and the imposed stress were then changed to 
T 2 and o 2 , respectively. After a steady state creep rate was attained, the 
test conditions reverted back to and T 1 if the specimen was still intact. 

Six step temperature cyclic tests were performed. As indicated in 
Table 3.3, two of the test specimens were subjected to fully reversed strain 
cycling with a strain range of ±0.4% at T]_. After cyclic saturation, the 
temperature was raised to T 2 and held at the peak strain for two hours. The 
test temperature was then lowered to T 1 , and the specimen was cycled again at 
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TABLE 3.3 


TEST MATRIX FOR THERMOMECHANICAL CONSTITUTIVE TESTING 

OF B1900+Hf 


(1.) Step Temperature Tensile Tests: 

Specimen Test # e, sec -1 T lf °C 

t 2 ,°c 

T 3 ,-c 

T,°C/sec 

K9 

1 

1 x 10“ 4 

538 

760 

538 

4.4 

K2 

2 

1 x 10 -4 

760 

982 

760 

4.4 

K1 

3 

1 x 10" 4 

982 

760 

982 

4.4 

K5 

4 

1 x 10" 4 

760 

982 

760 

4.4 

K4 

5 

4 x 10" 6 

760 

982 

760 

0.44 

K3,K11 

6 

1 x 10" 4 

760 

538 


4.4 

1A 

6A 

1 x 10 -4 

871 

649 


4.4 

(2.) 

Step Temperature Creep Tests: 




Specimen Test # 

o ^ 9 MPa 

T lf °C 

02*MPa 

t 2 ,°c 

T,°C/sec 

K8 

7 

690 

649 

427 

871 

4.4 

K6 

8 

690 

760 

215 

982 

4.4 

K7 

9 

215 

982 

690 

760 

-4.4 

(3.) 

Step Temperature Cyclic Tests: 




Specimen 

Test # , 

% Tj,°C 


t 2 ,°c 

R / e min\ 
K \ e max / 

T,°C/sec 

K17 

10 ±.4 

538 

hold at 

760 

-1 

4.4 



peak strain 




K16 

11 ±.4 

538 

hold at 

982 

-1 

4.4 



peak strain 




K12 

12 ±.4 

538 

±.4 

982 

-1 

4.4 


13 ±.4 

760 

±.4 

538 

-1 

4.4 

K13 

14 ±.4 

982 

±.4 

538 

-1 

- 4.4 
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TABLE 3.3 (CONTINUED) 


( 4 .) 

Specimen 

IB 

ID 

IE 

1C 


TEST MATRIX FOR THERMOMECHANICAL CONSTITUTIVE TESTING 

OF B1900+Hf 


Thermomechanical Cyclic Tests: 


Test # 

T,°C/sec 

T Range, °C 

e.sec * 

Ae ,% 

Loading Mode 

15 

2.8 

538-760 

lxlO' 4 

±.4 

in-phase 


1.85 

538-760 

lxlO" 4 

±.6 

in-phase 

16 

2.8 

538-760 

lxlO" 4 

±.4 

out-of-phase 


1.85 

538-760 

lxlO -4 

±.6 

out-of-phase 

17 

5.6 

538-982 

lxlO -4 

±.4 

out-of-phase 


3.7 

538-982 

lxlO' 4 

±.6 

out-of-phase 

18 

5.6 

538-982 

lxlO' 4 

±.4 

out-of-phase 

(Coated 

Specimens) 

3.7 

538-982 

lxlO' 4 

±.6 

out-of-phase 


47 


Temperature, T Strain 




FIGURE 3.6 STRAIN-TIME AND TEMPERATURE-TIME HISTORIES 
FOR STEP TEMPERATURE TENSILE TESTS. In 
some cases, more than one temperature 
change is planned. 
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FIGURE 3.7 STRESS-TIME AND TEMPERATURE-TIME HISTORIES FOR 
STEP TEMPERATURE CREEP TESTS 
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FIGURE 3.9 TMF LOADING HISTORIES OF TRIANGULAR WAVEFORM: (a) IN-PHASE LOADING, 

(b) OUT-OF-PHASE LOADING 







±0.4* strain until cyclic saturation was again obtained. In the other four 
cyclic tests, the specimens were cycled initially at ±0.4% strain at T 1 until 
cyclic saturation. The temperature was then changed to T 2 while cycling at 
±•4% strain continued. After cyclic saturation at T 2 occurred, the test tem- 
perature reverted back to Tj with continued cycling. 

Four thermomechanical fatigue tests were performed using transgranu- 
lar waveforms with the imposed temperature and strain being either in-phase 
[Figure 3.9(a)] or out-of-phase [Figure 3.9(b)]. In the former case, the tem- 
perature and strain imposed on the specimen increased and decreased linearly 
with time in similar manner reaching their respective maxima and minima at the 
same time. For the out-of-phase tests, temperature and strain varied linearly 
with time in opposite directions, e.g., a temperature increase was coupled to 
an imposed strain decrease (compression). The strain maxima occurred at the 
minimum temperature, while the maximum temperature occurred at the minimum 
strain. Both an in-phase and an out-of-phase cycle were performed for TMF 
cycling between 538 to 760°C. Only the out-of-phase cycle was performed for 
TMF cycling between 538 to 982°C. For examining coating effects, this out-of- 
phase cycle was also applied to the B1900+Hf specimen coated with yttrium- 
stabilized zirconia. All the four TMF tests were performed for strain ranges 
of ±0.4% and ±0.6% at a strain rate of 1x10"^ sec - *. 

3*4 Results: Experiment and Model 

3.4.1 Nonisothermal Tensile Tests 

Figure 3.10 shows the nonisothermal (step temperature) tensile tests 
conducted in the 760-982°C temperature range where thermal recovery is 
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important but without the possible complications of strain aging. The imposed 
strain rate was 1x10"^ sec - * and the heating or cooling rate was 4.4°C/sec. 
In one of the tests, the specimen (K3) was pulled isothermal ly at 760°C to 2 % 
strain. In a second test, the specimen (Kl) was pulled at 982°C to a strain 
of 1%, and the temperature was then decreased at 4.4°C/sec to 760°C. The tem- 
perature reached 760°C at a strain of 1.5%, and the specimen was then pulled 
isothermal ly at 760°C. The results in Figure 3.10 show that the nonisothermal 
specimen (Kl) exhibited similar stress-strain behavior at 760°C and appears to 
reach approximately the same value of saturation stress as the isothermal ly 
tested specimen (K3) indicating little thermal history effect. In a third 
test (Specimen K2), the test temperature was initially at 760°C and was then 
increased to 980°C after 1% strain. In this case, the stress-strain behavior 
of the nonisothermal specimen (K2) at 982°C was also consistent with that 
extrapolated from the isothermal curve for 982°C (Kl); i.e., the same value of 
stress at steady-state conditions. These results indicate that the steady- 
state tensile deformation behavior of B1900+Hf in the range of 760-982°C de- 
pends on the current test temperature but not on the thermal history. 

Before a comparison of model and experiment is presented for the re- 
sults shown in Figure 3.10, it is instructive to compare these nonisothermal 
test data, which are based on tubular specimens, with previously obtained iso- 
thermal data of round-bar specimens. Such a comparison, shown in Figure 3.11, 
reveals that the tubular specimens used for the nonisothermal tests are some- 
what weaker at 982°C but stronger at 760°C than the round-bar specimens. 
Figure 3.12 shows the nonisothermal stress-strain curves predicted by the 
Bodner-Partom model using material constants based on isothermal test data 
(round-bar specimens). Qualitatively, there is good agreement between theory 
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FIGURE 3.11 COMPARISON OF NONISOTHERMAL TENSILE RESULTS WITH ISOTHERMAL 
TENSILE DATA OBTAINED BASED ON ROUND-BAR SPECIMENS FROM 
BOTH THE BASE AND OPTIONAL PROGRAMS 
















MATERIAL ; B1900+MF : 

TEMPERATURE: 760-982°C : Bodner-Partora Model 
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FIGURE 3.12 BODNER-MRTOM MODEL PREDICTIONS OF NONISOTHERHAL TENSILE CURVES COMPARED TO ISOTHERMAL 



and experiment, as both the experimental and the predicted results indicate 
that absence of a thermal history effect on the steady-state or saturation 
stress at a given imposed strain rate. Quantitative comparisons of theory and 
experiment in Figures 3.13 and 3.14 indicate that the Bodner-Partom model 
overpredicts the stress at 982°C, underpredicts that at 760°C. The discrep- 
ancies appear because the model constants were determined from the isothermal 
baseline data which are slightly different than that of the tubular specimens 
used in the nonisothermal tensile experiments tested under the same 
conditions. 

Nonisothermal tensile tests were also conducted in the 538 to 760°C 
temperature range where strain aging effects could possibly occur [23]. In 
these tests, the specimens were strained at a constant rate of 10 sec at 
an initial temperature to a strain of - 2%. The test temperature was then 
changed to J 2 at a heating (or cooling) rate of 4.4°C/sec while constant 
straining continued. In all cases, the specimen temperature stabilized at 1 2 
at e> 2.5%. Figure 3.15 shows the nonisothermal tensile test results for 
B1900+Hf for the two cases: (1) straining initially at 538°C with the temper- 
ature subsequently raised to 760°C (Specimen K9), and (2) straining initially 
at 760°C with the temperature reduced to 538°C (Specimen K3 and Kll). The re- 
sults in Figure 3.15 indicate that the nonisothermal, plastic flow behavior of 
B1900+Hf at 760°C (Specimen K9) is in agreement with that extrapolated from 
the isothermal data (Specimen K3 and Kll). This behavior suggests that there 
is no thermal history effect on plastic flow, and it is consistent with the 
nonisothermal tensile data obtained in the temperature range of 760 to 
982°C. On the other hand, the plastic flow behavior of B1900+Hf at 538°C 
appears to depend, in part, on thermal history since different values of the 
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FIGURE 3.14 COMPARISON OF PREDICTED AND EXPERIMENTAL NONISOTHERMAL TENSILE CURVES 
FOR STRAINING AT e = lxlO -4 SEC" 1 WITH TEMPERATURE CHANGE FROM 982°C 
TO 760°C WITH ISOTHERMAL DATA 
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FIGURE 3.15. COMPARISON OF NONISOTHERMAL TENSILE RESULTS WITH MODEL PREDICTIONS FOR SPECIMENS 
TESTED IN THE 538-760°C TEMPERATURE RANGE 






saturation stress were obtained under isothermal and nonisothermal straining 
conditions. 

To further clarify the deformation behavior of Bl900+Hf in the 
temperature range of 538-760°C, an isothermal tensile test was performed at 
538°C under a strain rate of lxlO' 4 sec -1 (Specimen J3-1). Figure 3.15 shows 
that a large variation in the stress-strain behavior of B1900+Hf is observed 
between Specimen J3-1 and K9 at 538°C. The response curve for Specimen J3-1 
is in reasonable agreement with those of Specimens K3 and Kll for the portions 
at which T = 538°C. From Figure 3.15, it is clear that there is a large 
scatter in the isothermal data for B1900+Hf at 538°C. It is therefore not 
possible to conclude from these tests whether there is a thermal history 
effect associated with plastic straining of Bl900+Hf in the "strain aging" 
regime (538-760°C). This question will be raised again in relation to the 
nonisothermal cyclic tests. 

3.4.2 Nonisothermal Creep Tests 

The effects of thermal history on the steady-state creep behavior of 
B1900+Hf were investigated by performing step-temperature creep tests. In 
these tests, the specimens were subjected to a series of load and temperature 
histories designed to demonstrate whether or not the steady-state creep rate 
at a given stress and temperature depends on thermal history. Figure 3.16 
shows the nonisothermal creep results obtained by applying an initial stress 
of 215 MPa at 982°C (to Specimen K7) until a steady-state creep condition is 
reached (e = 1.3xl0' 6 sec -1 ). The temperature and the imposed stress were 
then changed to 760°C and 690 MPa, respectively, for which the steady-state 
creep rate became 7.5xl0' 7 sec' 1 . Subsequently, the stress was changed back 
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FIGURE 3.16 NONISOTHERMAL CONSTANT LOAD CREEP DATA OF B1900+Hf FOR TEST CONDITIONS 0 
o = 215 MPa AND T = 982°C WITH SUBSEQUENT STRESS AND TEMPERATURE CHANGES 
TO a - 690 MPa AND T = 760°C AND FINALLY TO a = 215 MPa AND T = 982°C 



to 215 MPa and the temperature to 982°C. The new steady-state creep rate then 
becomes 1.2xl0 -6 sec -1 , which is essentially the same as that obtained 
initially under similar test conditions. In Figure 3.17, the specimen (K6) 
was initially subjected to 690 MPa and 760°C which upon reaching steady state 
conditions, was changed to 215 MPa and 982°C. Examination of the test results 
in Figures 3.16 and 3.17 indicates that the steady-state creep rate of 
B1900+Hf depends on temperature and the applied stress, but not on the thermal 
history to within the experimental error. This conclusion also holds for 
creep behavior at 649°C. As shown in Figure 3.18, the creep rate for B1900+Hf 
at 690 MPa and 649°C, both before and after creep at 871°C and 427 MPa, was 
essentially zero. Because of the lack of thermal history effects, comparison 
of the noni sothermal creep data with those obtained isothermal ly in 
Figure 3.19 indicates good agreement between the two set of creep data, as well 
as with the Bodner-Partom and Walker models. 

3.4.3 Noni sothermal Cyclic Tests 

As in the nonisothermal tensile tests, the nonisothermal (step 
temperature) cyclic tests were performed at temperature ranges either in the 
thermal recovery regime (T > 760°C) or the "strain-aging" range (538 C < T < 
760°C). In Figure 3.20, the specimen was strained at ±.4% at 538°C until 
cyclic saturation. The test temperature was then increased to 760°C and was 
held at peak strain (.4%) for two hours. There was a slight stress relaxation 
during strain hold, as shown in Figure 3.21(a). After reducing the test 
temperature to 538°C, strain cycling was resumed at ±.4%, and the obtained 
saturated hysteresis stress-strain loops were essentially identical to those 
observed before the strain hold. The Bodner-Partom model predicted hysteresis 
loops, shown in Figure 3.21(b), which show fair agreement in the stress 
levels, both during cycling and strain hold. 
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FIGURE 3.17 NONISOTHERMAL CONSTANT LOAD CREEP DATA OF B1900+Hf FOR TEST CONDITIONS OF 
° = 690 MPa AND T = 760°C WITH SUBSEQUENT STRESS AND TEMPERATURE CHANGES 
TO a = 215 MPa AND T = 982°C 
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FIGURE 3.18 NONISOTHERMAL CONSTANT LOAD CREEP DATA OF B1900+Hf FOR TEST CONDITIONS OF 
o = 690 MPa AND T = 649°C WITH SUBSEQUENT STRESS AND TEMPERATURE CHANGES 
TO o = 427 MPa AND T - 871°C AND FINALLY TO a = 690 MPa AND T = 649°C 




FIGURE 3.19 COMPARISON OF NONISOTHERMAL CREEP DATA WITH 
ISOTHERMAL CREEP DATA AND MODEL PREDICTIONS 
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FIGURE 3.20 STRESS-STRAIN RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL 
CYCLIC SATURATION, STRAIN HOLD AT 760°C FOR TWO HOURS, AND THEN STRAIN CYCLING 
AT 538°C. The applied strain range is = ±.4% at e = 4x10"^ sec“l. 



Strain Hold at 760°C 
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FIGURE 3.21 STRESS-TIME RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL CYCLIC 
SATURATION, STRAIN HOLD AT 760°C FOR TWO HOURS, AND THEN STRAIN CYCLING AT 538°C 
The applied strain range is * ±.4% at £ = 4x10-4 sec"!; (a) experiment, and 
(b) model prediction. 
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FIGURE 3.21 (CONTINUED) STRESS-TIME RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL 

CYCLIC SATURATION, STRAIN HOLD AT 760°C FOR TWO HOURS, AND THEN STRAIN 
CYCLING AT 538°C. The applied strain range is - ±.4% at e = 4x10"^ sec*: 
(a) experiment, and (b) model prediction. 





A strain cycle similar to the one described in the previous 
paragraph was conducted on another specimen, but with strain hold at 982°C 
(Figure 3.22). Increasing the temperature for strain hold increased the 
amount of stress relaxation [Figure 3.23(a)]. Subsequent strain cycling at 
538 °C resulted in saturated hysteresis loops which were again essentially 
identical to those obtained before the strain hold. The Bodner-Partom model 
calculations for this strain cycle are shown in Figure 3.23(b). 

Additional cyclic tests were conducted under nonisothermal and con- 
tinuous straining conditions. Figure 3.24 shows the hysteresis loops for 
B1900+Hf strained at e = 4x10"^ sec"* and 538°C until cyclic saturation. The 
test temperature was then increased to 982°C, while the specimen was under 
concurrent strain cycling. After cyclic saturation occurred at 982°C, the 
test temperature was reduced to 538°C with the specimen under the same imposed 
strain cycling conditions. The subsequent cyclic saturated hysteresis loops 
at 538°C were found to be identical to those obtained before the nonisothermal 
histories were imposed (Figure 3.24). Rearranging the test temperature 
sequences to strain cycling at 982°C, 538°C, and then 982°C produced similar 
results with the cyclic saturated hysteresis loops for 982°C before and after 
thermal excursion to 538°C being identical (Figure 3.25). 

Figure 3.26(a) shows the nonisothermal cyclic results for the case 
in which the specimen (IF) was cycled at Ae = ±.4% initially at 760°C. After 
cyclic saturation, the test temperature was reduced to 538°C while the same 
strain cycling conditions were maintained. When cyclic saturation was 
achieved at 538°C, the temperature was then increased to 760°C, and strain 
cycling continued until saturation. The whole strain/temperature cycling 
histories were also repeated for Ae = ±.6%, and the results are shown in 
Figure 3.26(b). As shown in Figures 3.26(a) and (b), the cyclically saturated 
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FIGURE 3.22 STRESS-STRAIN RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL CYCLIC 
SATURATION, STRAIN HOLD AT 982°C FOR TWO HOURS, AND THEN STRAIN CYCLING AT 538°C. 

The applied strain range is = ±.4% at e = 4xl0"4 sec - *. 
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FIGURE 3.23 STRESS-TIME RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL CYCLIC 
SATURATION, STRAIN HOLD AT 982°C FOR TWO HOURS, AND THEN STRAIN CYCLING AT 538°C. 
The applied strain range is = ±.4% at e = 4x10"^ sec~l; (a) experiment, and 
(b) model prediction. 
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FIGURE 3.23 (CONTINUED) STRESS-TIME RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 538°C UNTIL 

CYCLIC SATURATION, STRAIN HOLD AT 982°C FOR TWO HOURS, AND THEN STRAIN 
CYCLING AT 538°C. The applied strain range is = ±.4% at e = 4xlO"^sec - l; 
(a) experiment, and (b) model prediction. 
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FIGURE 3.25 STRESS-STRAIN RESPONSE OF B1900+Hf SUBJECTED TO STRAIN CYCLING (1) AT S32°C UNTIL CYCLIC 
SATURATION, (2) AT 538°C UNTIL CYCLIC SATURATION, AND (3) AT 982°C UNTIL CYCLIC SATURA- 
TION. The applied strain range is = ±.4% at e = 4x10"^ sec - !. The hysteresis loops 
in (1) and (3) are essentially identical. 
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FIGURE 3.26 HYSTERESIS LOOPS OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 760°C WITH SUBSEQUENT 
TEMPERATURE CHANGES TO 538°C AND 760°C: (a) Ae = ±.4% and (b) Ac = ±.6%. The 

temperature changes were made after cyclic saturation was obtained 
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FIGURE 3.26 (CONTINUED) HYSTERESIS LOOPS OF B1900+Hf SUBJECTED TO STRAIN CYCLING AT 760°C WITH 

SUBSEQUENT TEMPERATURE CHANGES TO 538°C AND 760°C: (a) Ac = ±. 4 % and 

(b) Ac = ±. 6%. The temperature changes were made after cyclic satura- 
tion was obtained. 




hysteresis loops of B1900+Hf for T = 760°C do not depend on thermal history 
when strain cycled at Ae » ±.4% and Ae » ±.6% with temperature changes from 
760°C to 538°C and then back to 760°C. 

Comparisons of the nonisothermal cyclic data with previously 
obtained isothermal cyclic data shown in Figures 3.27, 3.28, and 3.29 for 
538°C, 760°C, and 982°C, respectively, indicate that the stress ranges of 
saturated hysteresis loops at a given strain range are independent of thermal 
history for specimens tested in either the strain aging regime (538°C) or in 
the regime where thermal recovery is important (982°C). These results are 
consistent with previous tensile and creep data obtained in the thermal 
recovery regime, which indicated that the saturation stresses and the steady 
state creep rates depend on the current temperature and loading conditions, 
but are independent of thermal history. 

3.4.4 Thermomechanical Fatigue Tests 

Two of the thermomechanical fatigue tests were performed in the 
temperature range of 538°C to 760°C. The strain rate was lxlO" 4 sec" 1 , and 
the imposed strain ranges were ±.4% and ±0.6% in each of the tests. In one 
test with Specimen IB, the imposed strain was in-phase with temperature as 
both strain and temperature increased or decreased linearly with time. For 
the other test, on Specimen ID, the imposed strain varied out-of-phase with 
temperature. Test results for both the in-phase and out-of-phase cycles in 
Figures 3.30-3.33 showed cyclic hardening at the hot end of the cycles (760°C) 
but not at the cold end (538°C). For the thermomechanical cycles at ±.6% 
strain range, the Bodner-Partom model predictions are in good agreement with 
experiment as shown in Figures 3.31 and 3.33. For cases at which Ae = ±0.4%, 
as shown in Figures 3.30 and 3.32, the Bodner-Partom model predicted 
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FIGURE 3.29 COMPARISON OF ISOTHERMAL AND NONISOTHERMAL CYCLIC 
DATA OF B1900+Hf AT 982°C 
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FIGURE 3.32 HYSTERESIS LOOPS OF B1900+Hf SUBJECTED TO OUT-OF-PHASE TMF 
LOADING UNDER A MECHANICAL STRAIN RANGE 0F« ±.37% IN THE 
538-760°C RANGE: (a) experiment, and (b) Bodner-Partom 

model calculation. 
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IMF loops are larger than that observed in the experiments. The discrepancies 
are primarily due to the small strain values (±.4%) in the near yield region 
of the stress-strain curves, which is very sensitive to specimen variation and 
small differences in the values of the model constants. An interesting 
computational result is that the constitutive model does indicate some degree 
of cyclic hardening at the hot end of the cycle. This asymmetrical cyclic 
hardening is less pronounced than the test results and saturates more rapidly 
at the higher strain range. 

The third thermomechanical fatigue test was an out-of-phase TMF 
cycle conducted in the temperature range of 538°C to 982°C at ±.0.4% and ±0.6% 
strain ranges. Under these conditions, thermal recovery is present at the hot 
end of the cycle (982°C), while strain aging might be present at the cold end 
(538-760°C) . The type of unusual, asymmetric cyclic hardening observed in the 
two previous TMF cycles was not present in these strain cycles as shown in 
Figures 3.34 and 3.35. Comparison of the predicted curves with experimental 
results indicates good overall agreement between theory and experiment with 
slight underprediction of the stress near the cold end of the cycle (538°C). 
It should be noted that all the TMF calculations shown in Figures 3.30-3.35 
were based on the Bodner-Partom model without any strain aging terms and which 
used material constants obtained from isothermal data. 

The fourth TMF test was conducted on a B1900+Hf specimen (Specimen 
1C) that was plasma sprayed with a yttrium-stabilized zirconia coating. The 
applied TMF cycle was identical to the out-of-phase cycle shown in 
Figures 3.34 and 3.35. The stress-strain response of the coated specimen, 
shown in Figures 3.36(a) and 3.36(b) for Ae = ±0.4% and ±0.6%, respectively. 
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FIGURE 3.34 HYSTERESIS LOOPS OF B1900+Hf SUBJECTED TO OUT-OF-PHASE TMF 
LOADING IN THE 538 TO 982°C RANGE AND A MECHANICAL STRAIN 
RANGE OF ±.37%: (a) experiment, and (b) Bodner-Partom model 

calculation . 
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FIGURE 3.35 HYSTERESIS LOOPS OF B1900+Hf SUBJECTED TO OUT-OF-PHASE TMF 
LOADING IN THE 538 TO 982°C RANGE AND A MECHANICAL STRAIN 
RANGE OF ±.58%: (a) experiment, and (b) Bodner-Partom 

model calculation. 
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Mechanical Strain (X) 

(b) 

FIGURE 3.36 EXPERIMENTAL HYSTERESIS LOOPS OF B1900+Hf WITH YTTRIUM-STABILIZED 
ZIRCONIA COATING SUBJECTED TO OUT-OF-PHASE TMF LOADING IN THE 
538 TO 982°C RANGE: (a) Ae = ±.4X, AND (b) Ae = ±.6X. 




is identical to that of the uncoated specimen (Specimen ID). In Figure 3.36, 
only the thickness of the metal substrate was used to compute the stress from 
the measured force. These results indicate that the ceramic coating did not 
support any loads under the imposed conditions. 

Several interesting observations were made with regard to the coated 
specimen: (1) cold spots were observed during TMF testing, and (2) bulges of 
the ceramic coating were observed after the TMF test. The cold spots and 
bulges could be the result of coating delamination along the coating/substrate 
interface. Microcracks were also observed on the ceramic coating, as 
illustrated in Figure 3.37. 

3.5 Microstructure, Deformation Modes and Dislocation Structure 

Transmission electron microscopy (TEM) was performed on selected 
B1900+Hf specimens to determine the influence of load and thermal history on 
the deformation structure of the alloy. The objectives of these efforts were 
as follow: (1) to identify the deformation mechanisms responsible for 
possible thermal -hi story effects in the B1900+Hf material, (2) to identify 
microstructural changes, if any, associated with particular thermomechanical 
paths, and (3) to provide information for relating the unified constitutive 
approach to the underlying deformation mechanisms and dislocation structures. 
As discussed earlier, thermal-history effects in the inelastic deformation of 
metals have generally been attributed to strain aging and other microstruc- 
tural influences. However, there is evidence in the literature which in- 
dicates that N i ^ A 1 exhibits anomalous deformation behavior in the intermediate 
temperature range (538-760°C) [16J. Since B1900+Hf contains approximately 60% 
(±5%) volume fraction of Ni^Al, the possibility exists that thermal-history 
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FIGURE 3.37 MICROCRACKS OBSERVED IN THE COATED B1900+Hf 
SPECIMEN SUBJECTED TO AN OUT-OF-PHASE TMF 
LOADING 
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effects in B1900+Hf could be caused by Ni 3 AI rather than by strain aging or 
changes in the microstructure. The TEM efforts were therefore intended to 
provide results for testing various proposals concerning thermal -hi story 
effects in the plastic deformation of metals. 

For the reason described above, the TEM efforts at SwRI were focused 
primarily on tests in which the specimens were deformed under nonisothermal 
conditions. In the companion Isotropic Life Prediction Contract, Moreno and 
his co-workers [24] at Pratt-Whitney Aircraft have examined the microstructure 
and the dislocation structures of the same B1900+Hf alloy deformed under 
isothermal conditions. In order to develop a more complete picture of the 
deformation behavior of B1900+Hf, important results from Moreno and associates 
for the isothermal B1900+Hf specimens will also be briefly reviewed. The 
results for the nonisothermal ly deformed specimens will then be presented and 
integrated with those of the isothermal specimens. 

Preparation of the TEM foils from the deformed B1900+Hf specimens 
involved sectioning and mechanical polishing of the specimens into thin 3 mm 
diameter disks followed by electropolishing of the disks into thin foils. 
Final thinning of the disks was performed in a Fischione jet electropolisher 
with a solution of 300 ml methanol, 220 ml ethylene glycol, and 30 ml 
perchloric acid maintained at -38°C. Following thinning, the precipitate and 
dislocation structure of the foils were characterized in a Philips 301 
transmission electron microscope. 

3.5.1 Microstructure 

As reported in a previous NASA report, the B1900+Hf material used in 
the present program was part of a single heat, designated W-0098, obtained 
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from Certified Alloy Products, Inc., Long Beach, California. The chemical 
composition, casting conditions, and heat-treatments of the B1900+Hf test bars 
were reported [2]. As indicated in the micrographs shown in Figure 3.38, the 
microstructure of B1900+Hf consists of mostly cuboidal y' in a y matrix. 
The Y ' size is approximately 0.6 ym in the as-cast material and is 0.9 ytn in 
the heat-treated condition. The volume fraction of y' is approximately 60%. 

Transmission electron microscopy performed at SwRI on deformed, 
nonisothermally tested B1900+Hf specimens has revealed three types of y' 
morphologies, including cuboidal y', spaghetti (elongated) y', and 
spherical y' of two different sizes which are shown in Figures 3.39(a), (b), 
and (c), respectively. These four types of y' were observed in all of the 
nonisothermally deformed B1900+Hf specimens. They were all sheared by 
dislocations with the highest dislocation density being in the y phase or at 
the y/y' interface. In most cases, the y' morphology, as revealed by TEM 
replicas shown in Figure 3.38, was of the cuboidal form. 


3.5.2 Dislocation Structures 
3. 5. 2.1 Isothermal Specimens 

In the Fatigue Life Prediction Contract, Moreno and his co-workers 
[24] examined the dislocation structure in B1900+Hf specimens deformed under 
isothermal loading. With the permission of the NASA manager. Dr. G. Halford, 
for the Fatigue Life Contracts, the results of Moreno, et al. [24] are 
summarized in this section for the purposes of comparing them with the 
nonisothermal results. For tensile specimens deformed at 538-1093°C at a 
strain of 1-1.5%, the dislocations in B1900+Hf were found mostly in 
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Fully Heat Treated y'= 0.9 n „ 10 f* m f 

FIGURE 3.38 GAMMA PRIME (y') SHOWS INCREASE WITH HEAT TREATMENT 
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FIGURE 3.39 THE MORPHOLOGIES OF y' OBSERVED IN DEFORMED 
B1900+Hf: (a) cuboidal y'» (b) spaghetti y'» 

(c) spheroidal y'. 
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the y phase and at the y h' interface forming a cellular network around the 
Y ' precipitates. Figure 3.40 shows that the dislocation network is very dis- 
tinctive in specimens deformed at temperatures above 760°C, but is less dis- 
tinctive at 538°C. In general, the dislocation density in the cellular net- 
work decreases with increasing test temperature which is probably due to 
static thermal recovery. Microtwins are also present in specimens tested at 
760°C. 

The dislocation structure found in isothermal creep specimens also 
consisted mostly of cellular networks surrounding the y' • This is illustrated 
in Figures 3. 41(a), (b), and (c) for creep specimens at (1) 982°C under a 
stress of 283 MPa for 4.1 hours, (2) at 871°C under a stress of 427 MPa for 
20.6 hours, and (3) at 760°C under a stress of 600 MPa for 134.8 hours, 
respectively. The creep strain for these specimens was approximately 3%, and 
the reduction in area was 2. 4-4. 9%. Note also that dislocation segments are 
present in the y' of these specimens. As in the tensile specimens, microtwins 
are observed in specimens subject to creep at 760°C. 

The dislocation structures of low-cycle fatigue specimens tested 
isothermal ly at 538 and at 871°C are shown in Figures 3.42 and 3.43, respec- 
tively. In both cases, the imposed strain ranges were 0.25%. The dislocation 
structures consist mostly of dislocation segments in y and at the y/y' inter- 
face. Dislocation segments in y' are also observed. Distinctive cellular 
dislocation networks, however, do not occur in these specimens primarily 
because the applied strain ranges are almost completely elastic. 

3. 5. 2. 2 Nonisothermal Specimens 

Six specimens of B1900+Hf were selected for dislocation structure 
characterization in the present program, and their loading histories are 
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A. 538°C, 1 .5% Strain 

B. 760°C, 1.3% Strain 
760°C, 4.5% Strain 

C. 871°C, 1.2% Strain 

D. 982°C, 1 .2% Strain 

E. 1093°C, 1.0% Strain 

TENSILE 



FIGURE 3 40 COMPOSITE TEM PHOTOS OF TENSILE TESTED B1900+Hf SPECIMENS 
AT VARIOUS TEMPERATURES, FROM MORENO, et al. [23] 
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FIGURE 3.41 DISLOCATION NETWORKS OBSERVED IN ISOTHERMAL CREEP 
SPECIMENS OF B1900+Hf: (a) creep tested at 982°C/ 

283 MPa for 4.1 hours; (b) creep tested at 871°C/ 
427 MPa for 20.6 hours; and (c) creep tested at 
760°C/600 MPa for 134.8 hours. Results are from 
Moreno, et al. [23] 


ORIGINAL PAGE IS 
OF POOR QUALITY 


98 


ORIGINAL PAGE IS 
OF POOR QUALITY 


□c 



99 


(a) creep tested at 982°C/283 MPa for 4.1 hours; (b) creep tested at 
871°C/427 MPa for 20.6 hours; and (c) creep tested at 760°C/600 MPa for 
134.8 hours. Results are from Moreno, et al. [23] 



FIGURE 3.42 DISLOCATION STRUCTURES AFTER BEING LCF TESTED AT 538°C 
(1000°F), 10 cpm, Ae = ±0.25% FOR 12,500 CYCLES, FROM 
MORENO, ET AL. [23] 
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FIGURE 3 43 DISLOCATION STRUCTURE AFTER BEING LCF TESTED AT 871°C (1600°F), 10 cpm, 
Ae = ±0.25% FOR 2,965 CYCLES, FROM MORENO, ET AL. [23] 


summarized in Table 3.4. One of these is an isothermal tensile specimen, 
three are nonisothermal tensile specimens, and the remaining two are 
thermomechanical fatigue specimens. These specimens were selected for the 
purpose of illustrating the effects of thermomechanical loading path on the 
dislocation structure. Figures 3.44-3.49 show the typical dislocation 
structures of the six specimens listed in Table 3.4. The main observations of 
the dislocation structures in the B1900+Hf specimens which are summarized in 
Table 3.5 are as follows: (1) there are extremely high dislocation densities 
in y and at the y/y' interface in specimens (Specimens J3-1, K3, K9, and 10) 
tested at the intermediate temperatures (538-760°C); (2) there are no large 
differences in the dislocation structures between the isothermal and noniso- 
thermal tensile specimens except that there are more straight dislocations in 
the y ' of Specimen J3-1 than in Specimen K3; (3) the dislocation density is 
reduced in specimens exposed to temperatures > 760°C (Specimens K2 and IE) 
because of thermal recovery; (4) straight dislocations are observed in y' in 
Specimen ID which was thermomechanical ly cycled between 538 and 760°C, while 
curved (wavy) dislocations are observed in Specimen IE which was thermomechan- 
ically cycled between 538 and 982°C. It is also important to note that all 
specimens exhibited extremely dense dislocations in the y and in the y/y' in- 
terface cellular networks surrounding the y" which is consistent with the pre- 
vious observations by Moreno, et al. Stacking faults were occasionally ob- 
served in B1900+Hf , but microtwins were not observed in any of the specimens 
tested. 

The Burgers vectors of the straight and curved dislocations in 
the y' of Specimens ID and IE were identified using the conventional two-beam 
extinction technique. In this technique, the two orientations of the g vector 
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TABLE 3.4 


B1900+Hf SPECIMENS SELECTED FOR MICROSTRUCTURAL CHARACTERIZATION 
USING TRANSMISSION ELECTRON MICROSCOPY 


Final Specimen 

Specimen Test Temp, °C Strain History Temp, °C 


J3-1 

538 

Tensile test 

538 

K2 

760/982/760 

Step temperature 
tensile test 

764 

K3 

760/538 

Step temperature 
tensile test 

538 

K9 

538/760 

Step temperature 
tensile test 

740 

ID 

538 - 760 

TMF 

649 

IE 

538 • 982 

TMF 

649 
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FIGURE 3.44 TEM MICROGRAPH OF SPECIMEN J3-1 SHOWING 
COMBINATION OF STRAIGHT AND WAVY DIS- 
LOCATION IN y' AND EXTREMELY HIGH DIS- 
LOCATION DENSITY IN y AND AT y/y' 
INTERFACE 



FIGURE 3.45 TEM MICROGRAPH OF SPECIMEN K2 SHOWING 
LOW DISLOCATION DENSITY IN y, HIGH 
DISLOCATION DENSITY AT y/y' INTERFACE 
AND CURVED (WAVY) DISLOCATIONS IN 
Y' 
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FIGURE 3.46 TEM MICROGRAPH OF SPECIMEN K3 SHOWING 
EXTREMELY HIGH DISLOCATION DENSITY IN 
Y AND AT y/y' INTERFACE, AND RELATIVELY 
LOW DENSITY OF STRAIGHT AND WAVY DIS- 
LOCATIONS IN Y 1 



FIGURE 3.47 TEM MICROGRAPH OF SPECIMEN K9 SHOWING 
SIMILAR DISLOCATION STRUCTURES AS IN 
SPECIMEN K3 
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FIGURE 3.48 TEM MICROGRAPH OF SPECIMEN ID SHOWING 
EXTREMELY HIGH DISLOCATION DENSITY IN 
Y AND AT y/y' INTERFACE, AND HIGH 
DENSITY OF STRAIGHT DISLOCATIONS IN 
Y' 



FIGURE 3.49 TEM MICROGRAPH OF SPECIMEN IE SHOWING 
HIGH DISLOCATION AT y/y' INTERFACE, 

LOW DISLOCATION DENSITY IN SOME y, AND 
RELATIVELY LOW DENSITY OF CURVED DIS- 
LOCATIONS IN y* 
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A SUKWARY OF DISLOCATION STRUCTURES IN B190Q+Hf SPECIMENS TESTED 
UNDER NON ISOTHERMAL OR IMF CONDITIONS 
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(the diffraction vector) for which the dislocations became extinct were 
located, and the Burgers vector b of the dislocations was obtained from the 
cross product of the two diffraction vectors based on the principle that g • 
b = 0. Figure 3.50(a) shows that the Burgers vector for the straight dis- 
locations in Specimen ID is of the <101> type. These are screw dislocations, 
i.e., the dislocation lines are aligned parallel to the Burgers vector. On 
the other hand, the curved dislocations exhibited in Specimen IE, shown in 
Figure 3.50(b), were found to be of mixed character having both edge and screw 
components since the dislocation lines are neither normal nor parallel to the 
Burgers vector. 

An extensive review of dislocation structures in Ni 3AI in either 
single crystal form or in y/y' alloys has been presented by Pope and Ezz [161. 
Comparing the present results with those summarized by Pope and Ezz [16], it 
appears that the straight dislocations observed in Specimen ID are probably 
sessile <101> (010) dislocations which are formed as the result of cross slip 
from the (111) to the (010) planes. Since these dislocations are immobile, 
they are barriers to subsequent dislocations either on the (111) planes or on 
the (010) planes. In principle, this hardening mechanism can be operative in 
tension, compression or cyclic loading of specimens at the intermediate 
temperature regime (538-760°C), but appears to be somewhat more pronounced 
under thermomechanical loading as reflected by the greater amount of straight 
dislocations observed in Specimen ID. The anomalous hardening observed in TMF 
cycling in the 538-760°C shown in Figures 3.31 and 3.33 for Specimens IB and 
ID is probably the result of the interactions of mobile dislocations with the 
sessile (010) dislocations. 
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(a) 



(b) 


FIGURE 3.50 DISLOCATIONS IN THE y 1 OF TMF SPECIMEN 
TESTED AT: (a) 538-760°C (Specimen ID), 

and (b) 538-982°C (Specimen IE). The 
straight dislocations in Specimen ID ex- 
hibit screw character, while those in 
Specimen IE show mixed character. 
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The above hypothesis is supported by experimental observations made 
on Specimen IE. In particular, it was observed that removal of the sessile, 
straight screw dislocations in y' by increasing the hot end of the TMF cycle 
to 982°C results in the elimination of the anomalous hardening. It is well- 
known that the rapid decrease in strength (critical resolved shear stress) of 
Ni 3AI and of y/y' alloys at temperatures above 871-982°C has generally been 
attributed to the activation of (010) slip (25]. At these elevated tempera- 
tures, the critical resolved shear stress for the (010) plane is lower and 
thermal recovery is more prominent than in the intermediate temperature regime 
(538-760°C). The higher test temperature in Specimen IE is believed to lower 
the critical resolved shear stress for (010) slip so that the dislocations on 
that plane are mobile and more favored those on the (111) plane. Any sessile 
<101> (010) dislocations that are produced at the cold end of the TMF cycle 
could become mobile at the hot end (982°C) of the TMF cycle thereby reducing 
the amount of straight screw dislocations that are observed. Without these 
sessile dislocations in y', there would be no additional barrier to dis- 
location motion and therefore no anomalous hardening. 

Intense shear bands were also observed in Specimens 10 and IE. As 
shown in Figure 3.51(a) and 3.51(b), these shear bands consist of narrow bands 
of material containing an extremely high density of dislocations which extend 
across both the y and y' phases. The slip bands were oriented on the (110] 
directions and were composed primarily of | [110] edge dislocations. Occa- 
sionally, slip bands were observed that contained dislocations of mixed or 
primarily screw character, as shown in Figure 3.51(c). Observed only in the 
TMF specimens which have been cycled at ±. 6 % for several cycles, these shear 
bands were probably formed as the result of the breakdown of the dislocation 
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FIGURE 3.51 INTENSE SHEAR BANDS OBSERVED IN SPECIMEN IE 
ARE COMPRISED PRIMARILY OF EDGE DISLOCATIONS 
(a and b), AND A SHEAR BAND IN SPECIMEN ID, 
WHICH IS COMPRISED OF MIXED DISLOCATIONS (c) 
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cellular structure due to strain cycling. They are likely to lead to fatigue 
crack formation, since they tend to concentrate the plastic strain. 

3.5.3 Relation of Microstructural Examinations and Unified 

Constitutive Theories 

One of the important observations of the dislocation arrangements in 
isothermal ly, as well as in nonisothermally deformed B1900+Hf specimens, is 
the formation of cellular dislocation networks around the y' precipitates. 
That fewer dislocations are observed within the y' can be rationalized on the 
basis that the phase ordered structure of y' requires dislocations within it 
to occur as pairs separated by an anti-phase boundary of relatively high 
energy [161. To minimize the total energy associated with the dislocation, 
dislocations which shear the y' are stored and rearranged in the y and at the 
y/y' interface into a cellular network. Such dislocation cell structures 
appear to form at relatively low strain levels (a few percent) and appear to 
be the stable configuration for specimens tested under different loading paths 
including isothermal tensile, isothermal creep, nonisothermal tensile, and TMF 
loading. The theoretical implication is that the saturation stress concept 
postulated in many unified constitutive models can be considered to be the 
flow stress associated with the saturated dislocation configuration, i.e., the 
stress required to move dislocations to and from the walls of the dislocation 
cellular network. Since the saturated dislocation structure is independent of 
the thermomechanical loading path, the saturation stress and the limiting 
values of the internal hardening variables should therefore be independent of 
the thermomechanical loading path. The values of the saturation stress 
should, however, depend on temperature and strain rate indicating the 
dependence of the critical resolved shear stress on these quantities. 
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The dislocation structure in specimens which have not reached the 
saturated state condition depends on their previous deformation history, and 
can be quite different in appearance from that of the saturated configuration. 
For example, the dislocation structures associated with strain cycling at low 
strain levels in Figures 3.42 and 3.43 consist mainly of dislocation segments 
in y' without a very distinctive cellular dislocation network in y and at 
the y/y' interface. It appears that the dislocation structure evolves with 
plastic deformation from a low density state to a saturated state with in- 
creasing inelastic strain or plastic work. The evolution equations for the 
hardening variable that are used in the unified approach are continuum re- 
presentations of the evolution of dislocation structures in the deformed 
metals. The rate of hardening corresponds on the rate of dislocation inter- 
action and rearrangement into cellular structures. For B1900+Hf and Mar-M247, 
the saturated configurations can be attained at a few percent strain. 

The rearrangement of dislocations into cellular structures im- 
plicitly suggests a dynamic recovery process since it is a means to reduce the 
interaction energy between dislocations. The observation that the cellular 
structure can be attained at low temperature without the benefit of thermal 
recovery clearly indicates the presence of a deformation-induced recovery pro- 
cess. The use of dynamic recovery terms in the evolution equations for the 
hardening variables, leading to saturation conditions, is thereby justified. 

The experimental evidence also shows that thermal recovery reduces 
the dislocation density in the cellular structure. The static thermal re- 
covery terms in the hardening equations are the continuum representations of 
this process. At elevated temperature, steady-state is controlled by the 
balance of the hardening and static thermal recovery terms. Consequently, the 
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high temperature steady-state structure can be different from the saturated 
structure. In B1900+Hf, this difference is manifested as a reduction in the 
dislocation density in the cellular structure. In the unified models, it is 
represented by the fact that the steady-state values of the hardening vari- 
ables (isotropic and directional hardening) are generally lower than their 
respective limiting values due to the presence of the static thermal recovery 
terms. 

The cell size of the dislocation configuration tends to decrease 
with stress in many metals. For B1900+Hf and possibly other strengthened 
alloys, the dislocation network is formed as the result of storing disloca- 
tions in y and at the y/y 1 interface. The size of the dislocation cell struc- 
ture appears to be controlled by the y' size and is relatively independent of 
the stress level. This effect can be seen by comparing the dislocation cell 
size associated with the cuboidal y' with those of the spaghetti y', with the 
latter being larger because of the larger spaghetti y 1 size. 

Thermal history effects in inelastic deformation are generally 
considered to be induced by strain aging or other microstructural influences 
[8,9,18]. For the B1900+Hf and Mar-M247 materials studied in the present 
program, there are no indications that strain aging plays an important role in 
the constitutive response of these alloys under thermomechanical loading. 
There are, however, some indications that interactions of dislocations in the 
ordered y' phase might be somewhat dependent on the thermomechanical loading 
path. This may be due to the forming of sessile (010) screw dislocation 
in y 1 which can occur only at the immediate temperature regime and is most 
pronounced under thermomechanical loading. Even under those conditions, this 
additional "anomalous" hardening effect is relatively minor compared to 
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conventional hardening due to y' and dislocation cell walls. Ignoring the 
effects of the sessile (010) dislocations in the Bodner-Partom model has not 
been found to introduce any significant errors in the predicted TMF calcula- 
tions. 

Based on these micromechanical considerations, it is fair to con- 
clude that the basic framework of the unified approach appears to be adequate 
for the class of materials studied in this program. In particular, the use of 
one kinetic equation for describing all aspects of inelastic flow, including 
plastic deformation, creep, and stress relaxation is fundamentally valid and 
supported by the dislocation structures observed on the deformed specimens. 
The general form of evolution equations for the hardening variables is also 
consistent with the various features of the dislocation configurations. For 
materials with stable microstructures, it does not appear that the hardening 
mechanisms under isothermal loading conditions are significantly different 
from those under nonisothermal loading conditions. It is therefore not 
necessary to introduce additional hardening terms into the evolution equations 
to account for thermal history effects other than those representing the vari- 
ation of material constants with temperature. 

3.6 Discussion of Thermal History Effects 

The conclusion obtained here that nonisothermal behavior can be 
predicted on the basis of isothermal data is different from those of other 
investigations [8-10,26]. It is therefore of interest to compare and discuss 
the present findings with those of Pelloux and Marchand [26] and of Robinson 
and associates [8-10]. In a recent report [26], Pelloux and co-workers 
examined the thermomechanical cyclic stress-strain response of B1900+Hf in the 
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temperature range from 400 to 925°C. By comparing the values of the stress 
range as a function of the strain amplitude for the TMF cycles with the 
isothermal strain cycles, they observed differences in the hardening behavior 
and concluded that it was difficult to predict the cyclic stress-strain 
behavior under TMF conditions from isothermal data. Their conclusion was 
directed mainly to the observation that the stress-strain responses of 
B1900+Hf under TMF cycling were different from those under isothermal 
conditions. Those differences arose because they compared the cyclic response 
characteristic of B1900+Hf subjected to very different thermomechanical 
loading histories. Since inelastic deformation is history dependent, it is 
therefore not surprising that a different response was observed for the 
different deformation paths; such differences of cyclic response were also 
observed in the present result. Our contention is that it is possible to 
obtain material constants in a unified model using only isothermal data and to 
apply the model for thermomechanical loading conditions. Without such a 
model, comparison of stress-strain responses of different thermomechanical 
histories is not meaningful. 

In a series of papers [8-10], Robinson and associates expressed the 
view that inelastic deformation under nonisothermal loading conditions cannot 
be predicted based only on isothermal data. In particular, they reported ex- 
perimental data of step temperature cyclic tests of Hastelloy-X, which showed 
that the stress response associated with strain cycling at 424°C after initial 
cycling at 649°C was different from the isothermal response at 424°C. On the 
other hand, the stress response associated with step temperature cycling from 
427”C to 649°C was in good agreement with the isothermal data. The latter be- 
havior was observed in the step temperature cyclic tests of the B1900+Hf 
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material in this program. The lack of agreement between the isothermal and 
noni sothermal data of Hastelloy-X subjected to strain cycling with temperature 
changes from 649°C to 424°C was attributed to dynamic strain aging [8-10], 
which was recently confirmed by Miner, et al [27], Since dynamic strain aging 
is not a factor in the B1900+Hf and MAR-M247 materials, the conclusions 
reached in the present program do not contradict those of Robinson. More 
recent work by Bartolotta and Ellis [28] on the Hastelloy X under isothermal 
and thermomechanical cycling suggests some correlation is possible between 
isothermal and noni sothermal cyclic response. Based on our experience and 
that of Robinson and co-workers, it is fair to conclude that isothermal data 
are sufficient for evaluating all the model constants in the unified models 
for materials which do not exhibit strong dynamic strain aging, plasticity- 
induced precipitation, or large microstructural changes in the temperature 
range of interest. The question of whether isothermal data is sufficient for 
evaluating all the model constants in a unified model for materials which 
exhibit concurrent plasticity-induced precipitation (e.g., 316 stainless 
steels and the Hastelloy-X alloy) remains open. Further work is required in 
this area. 
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4.0 APPLICATION OF CONSTITUTIVE MODEL TO AN ALTERNATE MATERIAL (TASK L) 


The application of the Bodner-Partom model to the alternate mate- 
rial, Mar-M247, was conducted in three phases. In Task LI, the minimum number 
of monotonic tensile, creep, and isothermal cyclic constitutive tests (both 
uniaxial and biaxial) was conducted to generate a data base for determining 
the model constants. Task LI was completed last year, while Tasks L2 and L3 
were completed this year. In Task L2, three cyclic thermomechanical tests 
were performed using the same test parameters that were employed in Task K for 
the base material. The applicability of the Bodner-Partom model was demon- 
strated for Mar-M247 in Task L3 by performing computations for the thermo- 
mechanical test cycles. 

For these tasks, a total of twenty (20) Mar-M247 specimens were 
fabricated. Fourteen (14) specimens were used for isothermal tension, creep, 
and cyclic strain tests; three (3) specimens were used for isothermal biaxial 
tests and three (3) specimens for the TMF tests. 

4.1 Experiments 

4-1.1 Isothermal Tensile. Creep, and Cyclic Testing 

Tensile tests with strain rate jumps were conducted for Mar-M247 
using the procedures described in [6]. The tensile tests involved extending 
the specimen at an initial strain rate, e.., to a prescribed 2 % strain level 
(corresponding to the saturation stress level). After reaching the prescribed 
strain, four changes of strain rate were imposed. The tensile test matrix for 
Mar-M247 is shown in Table 4.1. In addition, isothermal cyclic and creep 
tests were performed for Mar-M247. As indicated in the test matrix, the 
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TABLE 4.1 


ISOTHERMAL TENSILE, CYCLIC AND CREEP TEST 
MATRICES FOR Nar-M247 

Temperatures, °C 

23 648 760 671 982 1093 

Mar-M247 X X X X X X 

£.j = IxlO" 4 sec~^, s-j = IxlO" 6 sec”\ i 2 = IxlO* 5 sec~\ 

= lxl O' 4 sec~^ , = 1x10"^ sec"^ 


CYCLIC TEST MATRIX FOR Mar-M247 


Temperature, °i 

C e-j sec"^ 

Ae-| s ±.25% 

Ae 2 = ± .4% 

Ae 2 ~±.6% 

23 

5x1 0‘ 4 

X 

X 

X 

648 

5x1 0" 4 

X 

X 

X 

871 

5x1 0 -4 

X 

X 

X 

982 

5x1 O' 5 

V 
» \ 

V 

/\ 

X 

1093 

5x1 O' 5 

X 

X 

X 


CREEP TEST MATRIX 

FOR Mar- 

■M247 

Temperature, °C 

a 


871 

600 

MPa 

982 

300 

MPa 

1093 

150 

MPa 
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cyclic tests were conducted at three strain range levels (0.25%, 0.4% and 
0.6%) and at five temperatures (23, 648, 871, 982, and 1093°C). Only three 
creep tests were performed under conditions listed in Table 4.1. 

The biaxial deformation behavior of Mar-M247 was examined using 
three nonproportional strain cycles. The first strain cycle was conducted at 
648°C and contained the strain history sequence: (1) 90° out-of-phase com- 

bined tension-torsion, (2) in-phase combined tension-torsion, (3) fully 
reversed torsion, and (4) fully reversed tension. The strain sequence was 
conducted under strain-controlled condition and was repeated on the same 
specimen at two effective strain ranges (±.4% and ±.6%). The second nonpro- 
portional strain cycle was identical to the first one, but was performed at 
982°C. The third strain cycle was a combined tension-torsion cycle conducted 
under 90° out-of-phase condition with strain hold at each of the strain 
peaks. The test temperature was 982°C, and the imposed effective strain range 
was ±0.4%. 

4.1.2 Nonisothermal Tensile and TMF Testing 

Two nonisothermal (step temperature) tensile tests were performed on 
the Mar-M247 specimens. In one case, the specimen was pulled at 760°C at e = 
1x10“^ sec - ^ until a strain of « 1.0% was reached. The test temperature was 
then changed to 649°C. Straining at 649°C was continued until a strain of 
= 2% was reached. The test temperature was then changed back to 760°C. All 
the temperature change was done while straining at e = lxlO - ^ sec - ^. The 
other nonisothermal tensile test was performed in a similar manner, but the 
applied temperature sequence was 649°C and 760°C. An out-of-phase thermo- 
mechanical cycle identical to that of Specimen IE (see Table 3.3) was also 
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performed for Mar-M247. The temperature range for the TMF cycle was 538- 
982°C, and the imposed mechanical strain ranges were ±.4% and ±. 6 % at a rate 
of ±1x10“^ sec"*. 

4.2 Correlation of Model and Experiment 

Comparisons of calculations based on the Bodner-Partom model and ex- 
perimental results were presented for isothermal tensile, creep, and cyclic 
tests of Mar-M247 in the annual report [6] of last year (1986). For complete- 
ness, some of these comparisons are summarized here together with the 
nonisothermal test results obtained recently. The same set of material con- 
stants, which were determined from isothermal tensile and creep data shown in 
Table 4.2, was used for all the isothermal and nonisothermal calculations. 
The temperature dependence of material constants with temperature are shown in 
Table 4.3. The procedures for obtaining the material constants were described 
in Ref. [6]. Determination of the essential constants at the lower tempera- 
tures (23-648°C) was not possible since strain aging effects obscured the 
normal strain rate dependence of the flow stress. As a consequence, the value 
of obtained at the higher temperatures was assumed to be constant, and the 
associated value of n at the lower temperatures was obtained from the satur- 
ated flow stress. 

Correlation of model and experimental tensile stress-strain curves 
for six test temperatures at e = 1x10“^ sec - * is shown in Figure 4.1. The 
unusual temperature behavior at the lower temperature appears to be due to 
strain aging effects which were directly incorporated in the hardening 
constants. This procedure can be done for a single strain rate but does not 
demonstrate the influence of strain aging on the applicable strain rate 
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TABLE 4.2 


BODNER-PARTOM MODEL CONSTANTS FOR Mar-M247 

Temperature - Independent Constants 
m-j = . 1 MPa' 1 

a-j = 0.0 

Z 1 = 1.64xl0 4 MPa 

r l = r 2 = 3 
D q = IxlO 4 sec' 1 


Temperature - Dependent Constants 

Temperature, °C 


Constants 

23 

648 

760 

871 

982 

1093 

n 

.556 

.556 

.556 

.522 

.467 

.415 

m 2 (MPa' 1 ) 

1.35 

.47 

.47 

1.35 

1.35 

1.35 

Z Q (MPa) 

1 . 6xl0 4 

1.6xl0 4 

1 . 6xl0 4 

1 . 3xl0 4 

1.2xl0 4 

*3* 

O 

’>< 

in 

o 

r— 

Z 2 (MPa) 

1.6xl0 4 

1 . 6xl0 4 

1 . 6x1 0 4 

1.3xl0 4 

1 . 2xl0 4 

1.05xl0 4 

Z 3 (MPa) 

5xl0 3 

5xl0 3 

8.63xl0 3 

8.63xl0 3 

8. 63xl0 3 

8.63xl0 3 

A-|=A 2 (sec" 1 ) 

0 

0 

0 

2xl0' 2 

.25 

2 

E (MPa) 

1 .9x10° 

1 .61xl0 5 

1.71xl0 5 

1 .49xl0 5 

1 . 29xl0 5 

9 . 7x1 0 4 
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TABLE 4.3 


TEMPERATURE DEPENDENCE OF THE BODNER-PARTOM MODEL CONSTANTS 

FOR Mar-M247 


n = .556 


n = -.1575 + 


781.95 

T+273 


Z Q = Z 2 = 1.6 x 10 4 MPa 


Z Q = Z 2 = -6.5617 x 10 3 + 2 ’ 3 t +273 ~ MPa 


10 


A 1 = A 2 = 3.3064 x 10 exp 


3.214 x 10 
T+273 


T < 823°C 
823°C < T < 1093°C 
T < 760°C 
760°C < T < 1093°C 
] sec" 1 T < 1093°C 
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STRESS (MPA) STRESS (MPA) 



STRAIN C'O 
(a) 



STRAIN C*) 

(b) 


FIGURE 4.1 CORRELATION OF THE BODNER-PARTOM MODEL WITH EXPERIMENTAL TENSILE 
CURVES AT SIX TEMPERATURES: (a) EXPERIMENT, AND (b) MODEL 

CALCULATIONS. 
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sensitivity. As a consequence, the proper strain rate sensitivity in the 
presence of strain aging would not be indicated. Figures 4.2(a), (b), (c), 
(d), (e), and (f) show comparisons of model calculations and experimental re- 
sults for the stress-strain response under variable strain rate conditions at 
23, 648, 760, 871, 982, and 1093°C, respectively. The poor correlations of 
the predicted strain rate dependence of the flow stress with the test results 
at 23 and 648°C are due to strain aging effects which are currently not in- 
cluded in the model. Comparison of model and experimental creep curves under 
an applied stress of 60 MPa at 871°C is presented in Figure 4.3. For cyclic 
loading, the Bodner-Partom model predicts the stress range at a given imposed 
strain range reasonably well, but the calculated hysteresis loops are slightly 

overly squared, as demonstrated in Figure 4.4. 

Correlations of model and experiment for biaxial loading under 
nonproportional 90° out-of-phase tension-torsion strain cycling at 648°C are 
shown in Figure 4.5. The imposed effective strain range is ± 0.4%. As shown 
in Figure 4.5, the a-/3x loci at 648°C are predicted by the model, but 
individual hysteresis loops are larger than those experimentally observed. A 
comparison of model prediction and experimental data for 90° out-of-phase 
tension-torsion strain cycling with strain hold at the strain peaks is shown 
in Figure 4.6. In this case, the model also slightly overpredicts the 
tension/compression stress range and the width of the associated hysteresis 
loop. 

Figure 4.7 shows the noni sothermal tensile test result for the case 
in which the specimen was pulled initially at 760°C. The test temperature was 
then reduced to 649°C at a cooling rate of 4.4 °C/sec after the mechanical 
strain reached *.75%. The test temperature reached 649°C at a mechanical 
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STRAIN (',) 
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FIGURE 4.2 COMPARISONS OF MODEL CALCULATION AND EXPERIMENT FOR 
TENSILE TEST WITH STRAIN RATE CHANGES: (a) 23°C, 

(b) 648°C , and (c) 760°C. 
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STRAIN (5) 

(f) 


FIGURE 4.2 (Continued) COMPARISONS OF MODEL CALCULATION AND EXPERIMENT 

FOR TENSILE TEST WITH STRAIN RATE CHANGES: 

(d) 871°C, (e) 982°C, and (f) 1093°C. 







MATERIAL ; MAR-M247 
TEMPERATURE: 871 C 
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CURVE OF Mar-M247 AT 871 °C UNDER AN APPLIED STRESS OF 600 MPa 





STRESS (MPA) STRESS (MPA) 



-660 -.400 -.200 -.000 . 200 . 400 .500 000 

STRAIN (O 



-1 00 - 800 - S00 - 400 -.230 -.000 . 200 . 406 «00 800 1.000 

STRAIN CO 
(b) 

FIGURE 4.4 COMPARISON OF MODEL CALCULATION AND EXPERIMENTAL CYCLIC DATA AT 648°C. 

The cyclic test was conducted under fully reversed condition at strain 
ranges of ±.25%, ±.4%, and ±.6%: (a) experiment, and (b) Bodner- 

Partom Model . 
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EXPERIMENT 


Mar-M24 


BODNER-PARTOM MODEL 


vtar-«247 


-300 
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FIGURE 4.5 COMPARISON OF MODEL CALCULATION AND EXPERIMENT FOR Mar-M247 

UNDER 90° OUT-OF-PHASE STRAIN CYCLING AT 648°C. The effective 
strain and strain rate are +.4% and 4x10“^ sec - ^: (a ) o-t 

plot, (b) o-e plot, and (c) t-y plot. 
















MPa o, MPa 


EXPERIMENT BODNER-PARTOM MODEL 



-500 -300 -TOO 100 300 500 
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FIGURE 4.6 


COMPARISON OF MODEL CALCULATION AND EXPERIMENT FOR Mar-M247 
UNDER 90° OUT-OF-PHASE STRAIN CYCLING AT 982°C WITH STRAIN 
HULD AT EACH OF THE STRAIN PEAKS. The effective strain and 
strain rate are +_4xl0 - ^ sec - ': (a) a-x plot (b) a-e plot 
and (c) t-y plot. 
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FIGURE 4.7 COMPARISON OF THE BODNER-PARTOM MODEL CALCULATED CURVE WITH THE EXPERIMENTAL 
RESULTS OF A NONISOTHERMAL TENSILE TEST FOR Mar-M247 WITH TEMPERATURE CHANGES 
FROM 760°C TO 649°C AND THEN REVERTING TO 760°C 




strain of =1.0%, and straining continued until e - 1.85%. The test tempera- 
ture was then increased to and maintained at 760°C until the mechanical strain 
reached 3%. The imposed mechanical strain rate was lxlO" 4 sec" 1 throughout 
the test. Figure 4.7 also shows the Bodner-Partom model calculation obtained 
on the basis of material constants determined from isothermal data. Good 
agreement is observed between theory and experiment. In Figure 4.8, the test 
conditions were similar to the one in Figure 4.7, with the exception that the 
temperature sequence was 649°C and then 760°C. Good agreement between theory 
and experiment is also obtained in this case. Figures 4.9(a) and 4.10(a) show 
the experimental results for the out-of-phase TMF cycle for the mechanical 
strain ranges of ±.4% and ±.6%, respectively. The temperature range for the 
TMF cycle was 538-982°C, and the mechanical strain rate was lxlO" 4 sec -1 . The 
corresponding Bodner-Partom model calculations are shown in Figures 4.9(b) and 
4.10(b). The good agreement between theory and experiment indicates that the 
thermomechanical deformation behavior of Mar-M247 can be predicted on the 
basis of isothermal data only. 

4.3 Discussion 

The objective of Task L was to demonstrate the general applicability 
of the Bodner-Partom model and the associated procedures for determining 
material constant, which were developed for the material B1900+Hf, to the 
alternate material, Mar-M247. In order to evaluate critically the proposed 
methodology, only a limited amount of test data (mainly isothermal tensile and 
creep data) were used for obtaining the material constants. The Bodner-Partom 
model used in this exercise did not include a strain aging term, and all the 
model calculations including those of nonisothermal and TMF loading were based 
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FIGURE 4.8 COMPARISON OF THE BODNER-PARTOM MODEL PREDICTION WITH THE EXPERIMENTAL RESULTS 
OF A NONISOTHERMAL TENSILE TEST FOR Mar-M247 WITH TEMPERATURE CHANGE FROM 649°C 
TO 760°C 



STRAIN (MPa) STRESS (MPa) 



FIGURE 4.9 COMPARISON OF EXPERIMENTAL RESULT WITH BODNER-PARTOM MODEL 
CALCULATION FOR Mar-M247 SUBJECTED TO AN OUT-OF-PHASE TMF 
CYCLE OF Ac = ±.4%: (a) experiment, and (b) Bodner-Partom 

model calculation 
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FIGURE 4.10 COMPARISON OF EXPERIMENTAL RESULT WITH BODNER-PARTOM MODEL 
CALCULATION FOR Mar-M247 SUBJECTED TO AN OUT-OF-PHASE TMF 
CYCLE Ac = ±.6%: (a) experiment, and (b) Bodner-Partom 

model calculation. 
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on isothermal material constants. The good overall agreement between model 
and experiment for uniaxial and multiaxial strain cycling, and nonisothermal 
tensile and thermomechanical loading confirms the conclusion previously ob- 
tained for B1900+Hf : (1) material constants for the Bodner-Partom model can 
be determined based only on isothermal tensile and creep data; (2) no addi- 
tional terms are required in the evolution equations of the hardening vari- 
ables for representing thermal-history effects. The possible influence of 
thermal history on inelastic flow can be adequately obtained by integrating 
the constitutive equations for the imposed thermal and loading histories. 
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5.0 SUMMARY 


The concept of unified constitutive equations for elastic- 
viscoplastic material behavior has been in existence for twenty years, and the 
subject has undergone extensive development and application during the past 
ten years. This program has shown it to be a valid and reliable working tool 
for some of the more sophisticated engineering applications of materials where 
strain rate and temperature-dependent inelastic deformations are experienced. 
This conclusion is based upon: (1) the demonstrated capability of two avail- 
able unified constitutive models to correlate a wide ranging set of both uni- 
axial and multiaxial test data on nickel-base alloys at high temperature and 
(2) the efficient implementation and use of these constitutive equations in a 
general purpose finite element structural code. 

To summarize, the main accomplishments of this NASA-sponsored re- 
search program are: 


1. The unified approach for elastic-viscoplastic constitutive 
equations has been shown to be a valid and workable method for 
representing and predicting material response characteristics 
in the ranges where strain rate and temperature-dependent in- 
elastic deformations are experienced. 

2. A sophisticated test program on two high temperature alloys has 
been conducted over a wide temperature range for a variety of 
uniaxial and in- and out-of-phase multiaxial loading condi- 
tions, including thermomechanical cycling. The results of the 
test program provide a data base for evaluation of the repre- 
sentational and predictive capability of constitutive equa- 
tions. 

3. The program has motivated further development of the two sets 
of constitutive equations under examination. 

Walker Equations : modification of the governing kinetic equa- 

tion; some modifications to the hardening evolution equations; 
demonstration of temperature dependence of material constants. 
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Bodner-Partom Equations: development of procedure for deter- 

mining material constants from monotonic loading test data; 
development of method for including additional hardening due to 
nonproportional multi -axial loading; demonstration of tempera- 
ture dependence of material constants. 

4. Explicit procedures were developed for determining all material 
constants from isothermal test data. 

5. The two unified constitutive equations under examination have 
been incorporated into the MARC finite element code. A number 
of numerical exercises have been run including thermomechanical 
cycling histories on components of the hot section of gas 
turbine engines. 

6. Particular material effects have been examined in the test 
program with regard to their possible inclusion, if required, 
in the constitutive equations. These consist of: 


a. 

additional hardening due 
loading; 

to nonproportional 

multi axial 

b. 

thermal history effects; 



c. 

strain aging effects; and 



d. 

effects of thin coatings 
response. 

and grain size on 

mechanical 


7. Comparisons have been made of features of the two sets of con- 
stitutive equations with regard to alternative methods of 
describing directional hardening and measures of hardening. 
However, no decisive preferences can be stated. 
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